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Pesticides used improperly can be injurious to human beings, animals, 
and plants. Follow the directions and heed all precautions on labels. Store 
pesticides in original containers under lock and key—out of the reach of children 
and anirnals—and away from food and feed. 

Apply pesticides so that they do not endanger humans, livestock, crops, 
beneficial insects, fish, and wildlife. Do not apply pesticides where there is 
danger of drift when honey bees or other pollinating insects are visiting plants, 
or in ways that may contaminate water or leave illegal residues. 

Avoid prolonged inhalation of pesticide sprays or dusts; wear protective 
clothing and equipment, if specified on the label. 

If your hands become contaminated with a pesticide, do not eat or drink 
until you have washed. In case a pesticide is swallowed or gets in the eyes, 
follow the first aid treatment given on the label, and get prompt medical 
attention. If a pesticide is spilled on your skin or clothing, remove clothing 
immediately and wash skin thoroughly. 


NOTE: Some States have restrictions on the use of certain 
pesticides. Check your State and local regulations. Also, 
because registrations of pesticides are under constant 
review by the U.S Environmental Protection Agency, 
consult your local forest pathologist, county agnculture 
agent, or State extension specialist to be sure the intended 
use is still registered. 
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Foreword 


This report is published as a part of the USDA Forest 
Service program to improve aerial application of pesticides, 
specifically by using pesticides and delivery systems tailored 
to the forest environment. The program is conducted jointly 
by the Technology and Development Center, Missoula, MT, 
and the Forest Pest Management Staff, Washington Office 
at Davis, CA, under the sponsorship of State and Private 
Forestry. The MOD 6.0 version was cosponsored by the 
U.S. Army Dugway Proving Grounds and U.S. Army 
Chemical Research and Development Center, Research 
Directorate, Physics Division. 


Details of the aerial application improvement program are 
explained in two Forest Service reports, A Problem Analysis: 
Forest and Range Aerial Pesticide Application Technology 
(MTDC Rept. 7934 2804, July 1979, Missoula, MT) and 
Recommended Development Plan For An Aerial Spray 
Planning And Analysis System (Forest Pest Management 
Rpt. FFPM 82-2, February 1982, Davis, CA). 


A system of computer models has been developed to 
optimize spray program design and operation and assess 
environmental risk posed by aerial spray operations. 
AGDISP is one of these models, and this user’s manual 
covers the most recent version of AGDISP, MOD 6.0. 


AGDISP was developed under the leadership of the National 
Aeronautics and Space Administration with support by the 
Forest Service to include special forest spray conditions. 


When NASA’s Agricultural Aviation Program was terminated, 
the Forest Service continued to support and direct further 
development and improvements of the AGDISP model. 


The changes and additions that have been incorporated into 
MOD 6.0 are listed on page ii. 


The MOD 6.0 is operational on the USDA Forest Service 
Data General MV-15000 and is installed with GKS graphics, 
Digital Equipment BAX (with Tektronic 4025 graphics), and 
the IBM PC/XT/AT personal computers. 


Summary 


The AGDISP computer program predicts the motion of agricultural material 
released from aircraft, including the mean position of the material and the position variance 
about the mean as a result of turbulent fluctuations. Developed under sponsorship by 
NASA, U. S. Department of Agriculture Forest Service and U. S. Army Dugway Proving 
Grounds, this program operates efficiently and is user-friendly, with many of its features 
validated against wind tunnel and flight test data. This document is the Mod 6.0 User 
Manual to AGDISP, and includes a description of its current enhancements and operating 
instructions. This version of AGDISP is operational on the following computers: Data 
General (with GKS graphics); Digital Equipment Corporation VAX (with Tektronix 4025 
graphics); and IBM PC/XT/AT/compatibles (with proprietary plotting programs for dot 
matrix printers). 


Modification Level 


AGDISP program development has lead to the following released versions: 


Mod 2.0 Operational on a Control Data CYBER 175 at NASA-Langley. 


Mod 3.0 


Mod 4.0 


Includes: 


Graphics: 


All of the basic program development: fully rolled-up 
vortices or Betz roll-up, simplified terrain modeling, WAKE 
plot file entry, models for propeller, helicopter, crosswind, 
superequilibrium turbulence, canopy, vortex penetration into 
canopy, and material evaporation. Also includes a stand- 
alone program (AGLINE) to construct the equivalent 
Gaussian distribution. 


Tektronix 4025 and 401X terminals using NASA-Langley 
graphics software. 


Operational on a Univac 1108 at Ft. Collins (USDA Forest Service). 


Improvements: 


Additions: 


Graphics: 


Helicopter modeling with transition to rolled-up vortices; 
AGLINE calculations as a menu option in AGPLOT. 


Discrete crosswind velocity profile; nonzero deposition 
height; composite deposition plots (to 16 plot files); canopy 
penetration by helicopter downwash; plot option to plot 
material diameter time history. 


DISSPLA at Ft. Collins, with appropriate subroutine calls 
for Tektronix terminals. 


Operational on a VAX 11/785 at Dugway Proving Grounds (U. S. Army). 


Improvements: 


Additions: 


Graphics: 


Betz roll-up procedure and propeller model revised; 
equivalent Gaussian distribution selection criterion is a 
program decision involving material vertical velocity. 


Models for wide body effects, simple vortex circulation 
decay, jet engines, multiple powerplants and parameterized 
evaporation; default input file option; plot option to plot 
material vertical velocity time history. 


CALCOMP at Dugway, with appropriate subroutine calls for 
the pen plotter. 
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Mod 5.0 Operational on a Data General MV-15000 at Missoula, and on IBM PC/XT/AT 
personal computers and compatibles. 


Mod 6.0 


Improvements: 


Additions: 


Graphics: 


Revised solution procedure eliminates integration step size 
dependence on material decay time constant. 


Axial variation added to all models; ground sprayer; 
continuous deposition; canopy deposition; deposition on 
objects. 


GKS at Missoula; dot matrix printer output on IBM 
PC/XT/AT/compatibles. 


Additional features to enhance program usage. 


Improvements: 


Additions: 


Added inputs specifying up to 16 drop sizes in one run, and 
ability to compute both "to" and "fro" aircraft directions. 


Graphics extended to include flux through a vertical plane, 


volume and number median diameters, and coefficient of 
variation as a function of lane separation. 
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Modifications Since Mod 5.0 


The Program Card 0010 now contains only Maximum Time. Every AGDISP run 
is now considered a "full-plane" configuration. 


Last entries on cards 0025, 0029, 0055 and 0081 no longer need to be negative to 
signal last card. AGDISP now reads the input file once to establish the number of 
entries of each card type, then reads the input file a second time to acquire these 
values. 


AGDISP now internally proportions the aircraft weight on card 0023 to the biplane 
wings. Card 0023 should ALWAYS contain the assumed total aircraft weight. 


The jet engine exhaust plane, propeller blade plane, nozzle axial position, and wide 
body nose location are all measured relative to the trailing edge of the wing or the 
shaft centerline of the helicopter, now with the POSITIVE direction pointing 
TOWARD the nose of the aircraft. Thus, in most cases AGDISP expects the jet 
engine exhaust plane to be downstream of the wing (a negative number), the 
propeller blade plane to be upstream of the wing (a positive number), and the wide 
body nose to be upstream of the wing (a positive number). 


Card 0060 has been restructured. Now, the number of nozzles is entered as a 
single positive number; initial drop diameter is removed from this card and entered 
on new cards 0064; and the spray system flow rate is added to this card. Old 
AGDISP input files will NOT run with Mod 6.0 without revision. 


New card 0064 gives the initial drop diameter and mass fraction of each drop size to 
be solved (from one to 16). 


On card 0065, the entry of drop cutoff diameter DCUT has been changed to volatile 
fraction. 


New card 0090 has been added to permit a solution for the "fro" direction 
(AGDISP will always solve for the "to" direction). 


Gaussian ground deposition has been removed from the AGPLOT option list. 
Added options include computation of VMD and NMD, vertical flux through a 
specified plane, coefficient of variation and swath overlap pattern. 


Each plot generated in AGPLOT is automatically cycled into AGVIEW on personal 
computers, if the user is running with AGCTRL (preferred). 


Binary files are written in a slightly different format on personal computers, to 
reduce file size by ten percent. 


New card 0027 gives the Richardson number to correct the spread of the material 
from the nozzles for atmospheric stability. 
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List of Illustrations 


Schematic of the calculation plane of the AGDISP program. The 
aircraft is in level flight into the figure (the "to" direction). The 
positive x direction is out of the figure. 


Typical prediction and measurement of drop diameter variation 
with time as a consequence of evaporation (Dennison and 
Wedding 1984). 


Swirl velocity distribution of a rolled-up wake as a function of 
spanwise load distribution and as measured from the center of 
the resulting vortex. 


The composite velocity vector at an observation point found by 
summing the contributions of the aircraft vortex pair and its image 
pair below the surface. 


Schematic of the helicopter flow field model. 


Mean wind and root mean square turbulent fluctuations in a 
canopy (data from Wilson and Shaw 1977), with 

vx =k V(z,)/In(z,/Z,). 

Plant area fraction as inferred from the Chico almond orchard 
data in Newton, Barnes and Barry 1987. 

Inertial impaction efficiencies for targets discussed in Golovin 
and Putnam 1962. The top curve in each graph gives = 0; the 


middle curve gives @ = 100; and the bottom curve gives = 1000. 


Typical target normal vectors given as (x,y,z) components. 
AGDISP input file for Example Case 1. 


Mean trajectories for Example Case 1. The wing position is 
indicated by the double-wide solid line. 


Continuous ground deposition for Example Case 1. 
AGDISP input file for Example Case 2. 
Mean trajectories for Example Case 2. 


Mean trajectories (solid) and standard deviations around the 
mean trajectories (dashed) for Example Case 2. 
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Paths of the tip vortices (solid) and engine centerlines (dashed) 
for Example Case 2. 


Continuous ground deposition for Example Case 2. 
AGDISP input file for Example Case 3. 

Mean trajectories for Example Case 3. 

Continuous ground deposition for Example Case 3. 


Continuous ground deposition (nonvolatile) for Example Case 3. 
With the "normalized" scale selected, the area under the curve is 
also displayed. Card 0065 in Figure 6-9 shows a volatile fraction 
of 0.536, which means that ideally 0.464 of the released material 
will be nonvolatile. 


Volume median diameter at the ground for Example Case 3. The 
average drop size of 133.3 microns is also displayed. 


Number median diameter at the ground for Example Case 3. The 
average drop size of 132.9 microns is also displayed. 


Vertical flux through a plane located at y = -10 meters for Example 
Gases! 


Equivalent Gaussian distribution for Example Case 3. The 
corresponding figure of merit is also displayed. 


Crosswind velocity profile for Example Case 3. 


Diameter time history for material from the first nozzle in Example 
Case 3. Evaporation is nearly linear down to the nonvolatile 
diameter. 


U axial velocity time history for material from the first nozzle in 
Example Case 3. Material is released at the aircraft speed (49.6 
m/sec on card 0020 in Figure 6-9) and quickly adjusts to the 
background speed (from card 0028 in Figure 6-9) until the drop 
approaches the ground. 


V horizontal velocity time history for material from the first nozzle 
in Example Case 3. The drop is caught in the left tip vortex (as 
seen in Figure 6-10). 


W vertical velocity time history for material from the first nozzle in 
Example Case 3. The drop is caught in the left tip vortex (as seen 
in Figure 6-10). 


Deposition on the upstream side (-1,0,0) of a 0.075 meter diameter 


sphere placed in the Example Case 3 flow field at y = - 10 meters and 


z=2 meters. The target collection efficiency and final deposition 
level are also displayed. 
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Drift fraction for all material released in Example Case 3. The final 
fraction of material that remains drifting is also displayed. Figure 
6-9 (card 0065) indicates that a volatile fraction of 0.536 may 
evaporate. 


Drift fraction of nonvolatiles in Example Case 3. Figure 6-9 (card 
re, indicates that the nonvolatile fraction of released material is 
464. 


Coefficient of variation for the continuous ground deposition in 
Example Case 3. For a typical value of 0.3 (Teske, Twardus and 
Ekblad 1990), the ideal lane separation between aircraft flight lines 
is approximately 54 meters. 


Overlapped swath deposition pattern for Example Case 3, with a 
lane separation of 53.77 m as suggested in Figure 6-25. 


AGDISP input file for Example Case 4. 

Mean trajectories for Example Case 4. 

Continuous ground deposition for Example Case 4. 

Canopy deposition contours for Example Case 4. The near-zero 
contour is shown as solid, and the dashed contours are for 0.05, 
0.1, 0.15 and 0.2 gallons per acre. The maximum deposition is 
also displayed. 

Total canopy deposition for Example Case 4. 

Crosswind velocity profile for Example Case 4. The canopy 
extends to a height of 7.3 meters (from cards 0055 in Figure 
6-27); within the canopy the velocity is assumed linear. 


Plant area fraction profile for Example Case 4. 


Canopy deposition fraction for material from the third nozzle in 
Example Case 4. 


Drift fraction for Example Case 4. 


Program flowchart for operation of AGDISP and AGPLOT on 
personal computers. 


AGVIEW results with SAMPLE.PLT. The thick border defines 
the maximum plotting area. 


Turbulence factor in nonneutral atmospheric conditions referenced 


to neutral (Richardson number = 0; Bilanin, Teske and Hirsh 1978). 


AGDISP summary flowchart. 
AGPLOT summary flowchart. 
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Command Procedures for AGDISP on VAX 
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Nomenclature 


plant area fraction 

concentration 

effective canopy drag coefficient 
material or aircraft drag coefficient 
maximum body diameter 

material diameter 

target significant dimension 
inertial impaction efficiency 
gravity 

canopy height 

inertial impaction first parameter, Eq. (48) 
material mass 

ground deposition 

mean square turbulence level, equal to< uu >+<vv >+< ww> 
radius 

Reynolds number, Eq. (8) 
aircraft semispan 

wing planform area 

time 

aircraft thrust 

fluctuating local fluid velocity 


mean local fluid velocity 


aircraft flight speed 
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incremental axial velocity 

friction velocity 

fluctuating material velocity (u,v,w) 
mean material velocity (U,V,W) 
relative velocity, equal to| Uj - Vj| 
swirl velocity 

downwash velocity 

aircraft weight 

fluctuating material position (x,y,z) 
mean material position (X,Y,Z) 
canopy displacement thickness 
surface roughness 


reference height 

leaf capture efficiency 
circulation 

circulation at the wing centerline 
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1. Introduction 


The original development of AGDISP was motivated by a desire to determine how 
aircraft-unique wake and propulsion characteristics affect the ground deposition pattern of 
aerially released material. Material motion was computed by a Lagrangian formulation of 
the trajectory paths and the solution of ordinary differential equations. Equations for the 
ensemble averaged turbulent variance of the released material were developed using a 
locally isotropic spectral density function. Idealized flow field models were constructed for 
aircraft vortex wakes, helicopter flow fields, propeller swirl, terrain, crosswind, plant 
canopy and superequilibrium turbulence. The result of this work was the Mod 2.0 version 
of AGDISP as reported in Bilanin and Teske 1984 and Teske 1984. 


Enhancements to this program include additional models for wide body effects, jet 
engines and vortex circulation decay. Refinements to the program, notably to the helicopter 
model and graphical presentation of results, proceeded as the program was made 
operational at additional computer sites (Mod 3.0, Teske 1985; Mod 4.0, Teske 1986; Mod 
5.0, Teske 1988d, 1989b and 1990a). This document contains the latest model 
developments compiled in Mod 6.0, operational on the Data General computer at Missoula 
and IBM PC/XT/AT compatibles. 


This manual is organized as follows. Section 2 summarizes equation development 
and the models currently available in AGDISP. Section 3 discusses the operation of 
AGDISP and includes a listing of program restrictions, error and warning messages that 
may arise. Section 4 explains the input data options to AGDISP, and Section 5 explains 
the operation of the companion plotting package AGPLOT. Section 6 examines the test 
cases used to exercise AGDISP. Section 7 reviews the usage of AGDISP and AGPLOT 
on mainframe computer systems, while Section 8 reviews their usage on personal 
computers (including the use of AGEDIT, AGVIEW and AGCTRL). Section 9 highlights 
the selection of background turbulence level. The Appendices contain an explanation of 
each subroutine in the program, flowcharts of program structure and a listing of common 
variables. 
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2. Summary of Models 


A Lagrangian approach is used to develop the equations of motion of discrete 
material released from aircraft, with the resulting set of ordinary differential equations 
solved exactly from step to step. The formulation of the material equations, particularly the 
ensemble averaged approach used to develop turbulent properties of the material- 
atmospheric interaction, is detailed in Bilanin and Teske 1984. This section of the User 
Manual examines those variables and nomenclature that are a part of the AGDISP program. 


The aircraft is assumed to be in level flight near the surface, releasing material into a 
(y,z) plane normal to the flight direction, -x, as shown in Figure 2-1. AGDISP solves for 
the transport of this material within this plane, until it deposits upon the surface or is carried 
aloft by wind or vortical motions. Accuracy of predictions requires an accurate description 
of background environment (mean winds and turbulence). The need for accuracy must, 
however, be balanced by the requirement for computational speed and ease of calculation. 
Simulations need to be performed simply and repeatedly under a variety of background 
conditions, with minimum computer demands and setup requirements. When the need is 
justified, the solution for the background environment may be found by a larger, more 
detailed program, such as the UNIWAKE program discussed in Bilanin and Teske 1988 
and Teske 1988a. The mean velocity and turbulence profiles predicted by such programs 
may then be read by AGDISP to predict released material behavior. As an alternate 
approach AGDISP has been fitted with extensive simplified flow field options to permit the 
prediction of material motion in a wide range of idealized background conditions. These 
models are summarized in this section of the User Manual. 


The released material is assumed to be spherically shaped. The material flight path 
as a function of time after release is computed as the locations (X,Y,Z) for all material 
included in the simulation. The material velocity is denoted by (U,V,W). The material is 
affected by the background environment; however, the environment is assumed to be 
unaffected by the presence of the material. The interaction of the material with the 
turbulence in the environment creates turbulent correlation functions for the position and 
velocity, < yv > and < zw >; for the velocity variance, < vv > and < ww >; and for the 
position variance, < yy > and < zz >. The square root of these last two variables gives the 
horizontal and vertical standard deviations of the material motion about the mean described 
by Y and Z. 


MATERIAL EQUATIONS OF MOTION 


In AGDISP the released material is assumed to be influenced by two forces: weight 
and aerodynamic drag. The governing equations are derived in Bilanin and Teske 1984 
and repeated here for completeness: 


Mean Equations: 


d2X; 
dt2 
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Turbulent Correlations: 
G< xjxj > = 2 < xjvj > (3) 
G <xivj > =(< xij >-<xvj>)[1]+< vivj> (4) 
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In Eqs. (1) - (5), x V; and U; are the ensemble averaged ith components of 
material position, material velocity and local fluid velocity, respectively, while x;, v; and u; 


are the fluctuating ith components of material position, material velocity and local fluid 
velocity, respectively. 


Inherent in the material equations is a relaxation time T,, which is essentially the e- 


folding time for the released material to come up to speed with the local fluid velocity (for 
V; to approach and equal U;). The relaxation time is defined as 


Dp Pp 


Ae DP 
3 CE Vrel Pair (6) 


p 
where D,, is the diameter of the material, p,, is its density, p.;, is the density of air, V,,] is 
the relative velocity between the material and the local background, and Cp is the material 
drag coefficient. The semiempirical formula for Cp (Langmuir and Blodgett 1946) is 

Crp = 24( 1 + 0.197 Re0-63 + 0.00026 Re!-38 ) (7) 


where Re is the Reynolds number defined as 


Dy V. 
Reel rel : 
Vai (8) 
The material diameter Dp is affected by evaporation (Trayford and Welch 1977) as 
dis dDp fy 1 
Dp, dt eléti 
P 2 te | 1 oad (0) 


where T, is the e-folding evaporation time 


Dp 
(i 
B AO (10) 


and A@ is the wet bulb depression in deg C, with B defined as 
B = 84.76[ 1+ 0.27 Rel/2 | m2/sec - deg C : (11) 


In a recent series of tests (Dennison and Wedding 1984), the evaporation model 
described by Eq. (9) was compared with wind tunnel results. A typical comparison is 
shown in Figure 2-2, illustrating favorable agreement with the test data. 


Equations (1) - (5) cannot be solved without specifying relationships for the 
quantities < x;u; > and < ujv; >, the correlations of the material position and material 


velocity fluctuations, respectively, with the local fluid velocity fluctuation. These 
expressions are developed by integrating their ensemble averaged frequency spectra using a 
spectral density function for transverse velocity fluctuations in isotropic turbulence (von 
Karman and Howarth 1938). The results (detailed in Bilanin and Teske 1984) give 











2 t 
<xqu>= | 1) K+ 4 | (12) 
Peeee css 
eM Err (13) 
with 
T\2 T Ty \2 
-Cfh-2lef- 
eo i Loa Sails 
2 | sh 
imi 
Tr (14) 





where T, is the travel time of the material through a turbulent eddy of scale A, adjusted for 
the passive tracer limit 


eh 
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The quantity q? is the mean square turbulence level in the fluid. 


With the position and velocity information available for all released material at any 
time during the simulation, Eqs. (1) - (5) may be integrated exactly for the solution at the 
next time step (Teske 1988d). This step result serves as initial conditions for the next step 
integration, etc., through the entire simulation. The assumption that the background 
conditions Uj, < x;u; > and < u;v; > are constant across a time step permits solution of Eqs. 


(1) - (5) unrestricted by the size of the released material. Computational solution times are 


reduced significantly, and heretofore unsolvable material diameters (including the passive 
tracer limit, Dp = () are tracked by the equations. 


FLOW FIELD MODELING 


The behavior of the released material is intimately connected to the local background 
mean velocity Uj and turbulence field q2 through which the material is transported. 


AGDISP has been configured to accept user specified mean velocity and turbulence flow 
fields in neutral atmospheres, but also contains a number of simplified models for the flow 
field velocity and turbulence levels behind aircraft. These models are summarized below. 


Fixed-Wing Rolled Up Tip Vortices 


When an aircraft flies at a constant altitude, the aerodynamic lift generated by the 
lifting surfaces of the aircraft equals the aircraft weight. The majority of the lift is carried 
by the wings, and generates one or more pairs of swirling masses of air (called vortices) 
downstream of the aircraft. If the roll-up of this trailing vorticity can be approximated as 
occurring immediately downstream of the wing, then the mean velocity field that results 


may be simply characterized by the aircraft semispan s, circulation I and load distribution 
(Figure 2-3). 


Fixed-Wing Betz Roll Up 


When the wing loading cannot be approximated by a single vortex pair, or the 
aircraft is flying sufficiently fast that roll-up cannot be assumed to occur immediately 
downstream of the wing, then the Betz methodology (Betz 1933) may be employed. The 
solution procedure (documented in Bilanin and Donaldson 1975) relates the swirling 
velocity distribution in the vortices to the details of the wing spanwise load distribution by 
assuming that angular momentum is approximately conserved. For complicated wing 
planform shapes, a vortex lattice analysis must be invoked to generate the wing spanwise 
load distribution (Margason and Lamar 1971). In a simulation involving the Betz 
methodology, the vortex sheet rolls into the vortex as a function of time. Thus, the 
unrolled sheet will contribute to the ambient fluid velocity during the roll-up process. 


Vortex Circulation Decay 


The flow field generated by a simply loaded wing includes four vortices: the two 
from either wing tip, and image vortices below the surface (Figure 2-4). The local ambient 
velocity is the vector sum of the four swirl velocity components, each written as 


i 


i lea (16) 
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with its direction measured perpendicular to a line connecting each vortex centroid and the 
observation point. The observation point might be a material position or the location of 
another vortex (which also moves because of the other vortices present in the flow). An 


inviscid vortex of constant positive circulation strength I, descends toward the surface 
while increasing the separation distance between its companion vortex. 


In the atmosphere, turbulence acts to decay the vortex strength. A simple decay 
model, developed in Donaldson and Bilanin 1975, results in a vortex decay of the form 


P= exp (- 22] (17) 


S 


For vortices out of ground effect, a typical value for B = 0.41; while a detailed examination 
of Program WIND data in ground effect (Teske 1988b and Bilanin et al. 1989b) gives Bq = 
0.56 m/sec. 


Helicopter In Forward Flight 


The helicopter model includes both hover downwash and the tip vortex pair by 
partitioning the helicopter weight between the two effects as a function of time. The hover 
downwash model is taken from actuator disk theory for a propeller and may be written as 


F Wy = 2% Pair R2 wa (18) 
while the strength of the vortex pair may be found from 


where W, is the helicopter weight, wg is the downwash at the rotor plane, R is the radius 


of the rotor and I is the tip circulation strength. The function F is determined by the 
expression 


(20) 


where | is the forward advance ratio, © is the solidity and I’, is the solution to Eq. (19) for 


F=0. The constant k relates helicopter roll-up (around a circumference at the blade tip) to 

fixed-wing roll-up. Predicted results from a helicopter wake model (Bliss, Teske and 

Quackenbush 1984) suggest that the downwash flow field quickly transitions (within two 

blade revolutions) into a vortex pair. A typical value for k = 7. Comparisons of Eqs. (18) 

; neae aed detailed helicopter wake predictions give surprisingly good results (Teske 
a). 


The tip vortices influence the flow field everywhere, while the downwash affects 
only the spreading region beneath the rotor and the fluid column directly above the rotor. 
Within the boundaries of the rotor blade and the dividing streamline (Figure 2-5), the 


velocity is assumed to follow a stagnation point flow, while the vortex pair created by the 
helicopter is assumed to follow the dividing streamline. 


Material released ahead of the helicopter (spray boom forward) is assumed to 
encounter a streamline pattern similar to flow around a circular cylinder (the actuator disk 
downwash, Kuethe and Schetzer 1964) 


2 2R2y2 
is eg poet Ped Seger tie 
eeu sn a4 (21) 
2U,,R2xy 
Mee eA (22) 


until the material crosses the plane of the helicopter shaft centerline. 
Jet Engine 


The jet engine exhaust is modeled as a turbulent circular jet using the similarity 
analysis reported in Schlicting 1968. The results give the axial and radial velocities as 
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where T is the jet thrust and x is the downstream distance (corrected for the virtual origin). 
The radial velocity V,agiaj becomes part of the local fluid velocity, while the axial velocity 


eal is also used to evaluate the turbulence level 


: 2 
q2 = 0.2034 u2..1 (27) 
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where the proportionality constant is determined from the centerline decay of turbulence in 
a free jet (Wygnanski and Fiedler 1969). 


Propeller 


The propeller is modeled as an actuator disk, where the incremental velocity AU 
over the flight speed U,, is related to the actual thrust produced by the propeller 


T=2m pair R? AU(U,, + AU | (28) 


where R is the radius of the propeller. In steady flight the thrust equals the aircraft drag, so 
that 


Sele ee 


where S is the wing planform area and Cp is the aircraft drag coefficient. Combining Eqs. 
(28) and (29) to eliminate thrust, we obtain 





AU—1(—14 De 
com ee mt R2 





(ew 


(30) 


An alternate expression for €, defined in Eq. (26), may be obtained from Schlicting 1968 
as 


€ = 0.0256 R AU (31) 


Combining this equation with Eq. (26) gives the effective thrust level for the propeller as 
Tere = 2.5275 pair R% AUZ (32) 


to permit the use of Eqs. (23) - (27) for the propeller when Eq. (30) is substituted. The 
swirl velocity generated by the propeller is assumed to be linear in R out to its maximum 
value, and then zero for larger values of r. The resulting integration yields 


U> Cps 
V ip =a ee ee 
i rn QR3(U,, +AU | (33) 


where 71] is the propeller efficiency and Q is the propeller rotation rate. 


2] 


Wide Body Effects 


The presence of a significant fuselage shape may generate some additional mean 
velocity and turbulence in the vicinity of the aircraft. For simplicity, we assume that the 
fuselage cross-sectional area can be modeled as an axisymmetric body shape B(x), where x 
is measured from the nose to the tail. Slender body theory (Ashley and Landahl 1965) then 
gives the additional radial velocity contribution from the changing area of 


U 


co 





. = B’ 
Vradial 5 Hee (x) (34) 


where the prime denotes a derivative in x. 


The turbulence level, consistent with the axisymmetric wake created by the passage 
of the body, is found by matching the data in Chevray 1968 with a turbulence model 
(Lewellen 1977) and assuming a simple exponential decay 





pie ipl oN ty ees 

t= Ibody exp | 2(Rip | (35) 
where 

Dear. 2{x/D \1.43 

Abody = 0-003 Us, ( ti (36) 
and 

Ri2=D{ =LR\p 9: (37) 


is the wake radius to half-width, D is equal to the maximum diameter of the body and x is 
corrected for the virtual origin. 


Mean Crosswind 


In a neutral atmospheric surface layer the horizontal velocity follows a logarithmic 
profile 


- Tn z4)2 75) 
V(z) = V(zq) hears) (38) 


where V(z,) is the known horizontal velocity at a given altitude z,, and Z, is the surface 
roughness. Surface roughness is typically taken to be 1/30th of the actual physical 
roughness height. With a linear integral scale of turbulence (A = 0.65z), the turbulence 
level becomes 





V(x) i 

2 = 0.845 | —— _ 

4 In( z;/Zg ) (39) 
Superequilibrium Turbulence 


The detailed effects of turbulence are obtained by invoking superequilibrium 
turbulent transport theory (Donaldson 1973). Superequilibrium refers to the second-order 
closure turbulence transport model limit where the velocity correlations are able to track 
their equilibrium values. The nonlinear equations solved for the turbulence may be written 
in index form as 


ij 3 
i , : : q Pate ns 1M 4 39) SU ay 
[ Bim <untj> + Bj < vita > ]+ | < uu > 3 Oj 4 rat ary 
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Since the local mean flow gradients are known, the system of equations represented by Eq. 


(40) may be solved exactly for < uju; > to determine q~. 


Plant Canopy 


The canopy flow field model is an approximation of the second-order closure 
turbulence model discussed in Wilson and Shaw 1977. For a canopy of height h, and 


plant area fraction A(z), the mean crosswind velocity above the canopy is 


In(z/ Zo) 
V(z) = V(z-) —\— 
where z, is the displacement thickness of the canopy. Within the canopy the mean velocity 


and root mean square turbulence levels are assumed linear with height, in good agreement 
with data (Figure 2-6). The canopy alters the downwash stagnation flow below the 
helicopter by moving the effective surface up to z,. Wortices entering the canopy will have 


decaying circulation strength of the form 





ol ich 
To Vo 
1+ nae f(t) (42) 


similar to the atmospheric decay effect in Eq. (17). Here the function f(t) integrates the 
cumulative effects of A(z) on the circulation of the vortex 


prods 
f= af A(z) fa dz 
Ah JO dhe - Ah (43) 
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where Cc is the canopy drag coefficient, Ah is the penetration depth of the vortex and f Ais 
the fraction of the vortex within the canopy. 


Terrain 


Surface slope is modeled by assuming that the ground plane remains straight but 
may incline left or right. Any vortices in the flow field will have a modified image vortex 
system to maintain zero flow through the tilted surface. The helicopter downwash model 
and the position of the dividing streamline are also altered. Crosswind and canopy effects 
are assumed to remain parallel to the tilted surface. 


Ground Sprayer 


A ground sprayer is modeled by eliminating the presence of the aircraft vortices and 
being very specific about the direction and initial velocity of each nozzle. 


DEPOSITION MODELING 


As material approaches the surface, deposition begins, and continues until all 
material is deposited (if evaporation occurs, some of the material will be left in the 
atmosphere to drift). Ground deposition is computed by assuming that the concentration of 
material around the mean may be taken as Gaussian 


peewee ((z-ZP 
202 202 


c=—l exp 
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exp (44) 


where the released material is at position (Y,Z). 
Gaussian Deposition 


The material is assumed to deposit entirely at the point of surface impact. Here, 
Eq. (44) is integrated across all z values to give 


1 (y-Y} 
M = exp| - ————— 
plane Sari P| Bre 








(45) 


Continuous Deposition 


The material is assumed to deposit incrementally as it approaches the surface, and to 
continue depositing at its approach rate after impact. Here, Eq. (44) is integrated from far 
below the surface to the material location to give 


M ee ee 
plane ONDE p 








ee erfc ( A | (46) 


where erfc is the complementary error function. Ground deposition is obtained by 
summing all incremental contributions to Mplane as the integration proceeds. It may be 


seen that for material falling vertically toward the surface, the ground deposition pattern 
generated by Eg. (46) will be identical to the Gaussian deposition of Eq. (45). The use of 
Eq. (46) leads to more realistic ground deposition distributions. 


Canopy Deposition 


When released material traverses a canopy, deposition on the canopy reduces the 
amount of material available for deposition on lower levels of the canopy or on the surface. 
The amount of material captured on the canopy vegetation will be directly related to the 


capture efficiency f and the distance traveled by the material through the canopy. Thus, at 
any time increment At , the incrementally deposited material may be found from 


AC = B C A(z) Ad (47) 


where Ad is the distance traveled in At. The quantity AC reduces the remaining material C 
available for further deposition. 


Equation (47) contains A(z), the plant area fraction. This is essentially a discrete 
function of height through the canopy, and is the ratio of leaf area (space "occupied" by the 
tree) divided by the surface area allocated to each typical tree. Newton, Barnes and Barry 
1987 measured the leaf area for typical almond trees in the Chico orchard. Since the trees 
were spaced 8.2 meters apart, the surface area for each tree becomes 68 sq. meters. The 
resulting plant area fraction, discrete at seven heights through the canopy, is shown in 
Figure 2-7. 


Object Deposition 


Equation (44) may also be interpreted to recover the deposition on collection 
devices (cards, cylinders or spheres) placed at specific locations in the wake flow field. 
The total amount of deposition, as computed by Eq. (44), is reduced by the efficiency of 
the collection device, using the formulas and empirical data found in Golovin and Putnam 
1962. The impaction efficiency E (as plotted in Figure 2-8 for the three collection devices) 
depends on two parameters 


ae Pp Dé U 
9 Pair Vair Dt (48) 

6 Rs 9 Pair D; U 
Vair Pp (49) 


The significant target dimension D, is a typical cross-sectional distance used to characterize 


the size of the target (for a card, it is the smaller length; for a cylinder, it is the smaller of its 
diameter or length; and for a sphere, it is its diameter). Figure 2-8 demonstrates that as the 


target dimension increases, K decreases from Eq. (48), @ increases from Eq. (49), and the 
collection efficiency decreases (it is easier for a drop to get out of the way of a big target 
that significantly disturbs the flow field than a small target that does not). 


Material impaction also requires a knowledge of the orientation of the target relative 
to the solution coordinate system. For the purposes of this analysis, a normal vector must 
be supplied. This (x,y,z) vector defines the direction of a perpendicular line away from the 
target (in the case of a cylinder or a sphere this information isolates which half of the target 
is to be investigated). Some typical normal vectors are given in Figure 2-9. 


OVERALL MODEL VALIDATION 


The first significant comparison of AGDISP predictions with data was reported in 
Morris et al. 1982, where the influence of rolled-up tip vortices, propeller and crosswind in 
the mean and variance equations were favorably compared with a series of fixed-wing 
experiments. Later detailed comparisons (Teske 1988c, 1989a, 1989c and Bilanin et al. 
1987) demonstrate the favorable performance of AGDISP when compared with detailed 
field test data for both fixed-wing and helicopter experiments over a wide material size 
range. AGDISP models are also summarized in Bilanin et al. 1989a, and used to predict 
lane separation for Gypsy moth spraying in Teske, Twardus and Ekblad 1990. 
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Figure 2-3 Swirl velocity distribution of a rolled-up wake as a function of spanwise load 
distribution and as measured from the center of the resulting vortex. 
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Figure 2-5 Schematic of the helicopter flow field model. 


(92/72) uy / CZ)AX = *A 


yim “(2/61 MeYS pur UOSTIAA Woy eIep) Adouvo & UI SUOnENIONT JUsTNGImM sends uUvdUT JOOI PUL PUIM URI, 9-7 DINIL 





*A/  <MM>‘*A/ <AA> BAVA 
a ); (a2) 
OZ 0"! O 8 9 p 2 0 
O 0 
20 20 
<MM> 
acd | > 2 dbus = 
of = OO uv BlIfOud 
elu pawnssy— 
oa 9°0 90 
ne, ni 
QO Vale oy Oy 
go 4 go 2 
Go ae 
. a }YB198H be 
ri a kdouv) 
Z| cel 
® 18 
QO { I 
© 18 
bi | pI 


2-18 


"L861 Alreg pue soureg ‘uojMONY Ul Byep PIeYoIO puoUTTe OdIYD oy} WIOIY PoLlayur se uOMOVY wae jue[g /-7 N31 


Gaale, 
Se@ eee ~=—S STC 


Q 


CAD Lao Sri 





ae 


3341 GNOW1V 





2-19 


1.@ CARD 
0.8 
8-6 
0.4 
8-2 


“| Q 1 e 
LOGCK) 


Pom SNDER 





LOGK} 


1.9 SPHERE 


8.8 
8.6 
0.4 
B.2 
ape 
-1 1 2 
LOG) 


Figure 2-8 Inertial impaction efficiencies for targets discussed in Golovin and Putnam 
1962. The top curve in each graph gives @ = 0; the middle curve gives 6 = 


100; and the bottom curve gives 6 = 1000. 
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Figure 2-9 Typical target normal vectors given as (x,y,z) Components. 


3. Code Operation and Error Messages 


AGDISP really consists of two programs: AGDISP for establishing the desired 
background fields and computing the material trajectories, and AGPLOT for plotting the 
resulting solutions. These programs are configured as interactive programs reading their 
needed input data from disk files and writing output to the terminal. AGDISP reads an 
input file of data cards (valid cards are discussed in Section 4), displays information to the 
screen as the run proceeds, generates printed output in a separate file, stores trajectory 
information in plot files for subsequent plotting, and reads WAKE plot file results from a 
separate data file (if invoked). AGPLOT reads the desired material plot files and 
interactively questions the user to supply necessary scaling data for plot construction. 


AGDISP reads and processes the input data file, forcing termination of the program 
if errors are discovered or needed inputs are missing. Initialization follows: material 
locations and velocities, vortex strength, Betz roll-up, plant canopy, crosswind velocity 
profile, WAKE plot file. The Lagrangian equations are then integrated repeatedly until one 
of three termination conditions occurs: 

1. Maximum simulation time is reached; 

2. All of the released material has deposited on the surface; or 


3. All of the material has reached the evaporation cutoff diameter and are removed 
from the calculation (if invoked). 


Each integration step in turn includes the following: 
1. Equation system solution of the equations of motion; 


2. Updating the background mean velocity and turbulence fields by the WAKE 
plot file, Betz roll-up, vortex motion, or plant canopy; and 


3. Incremental saving of the results for subsequent plotting by AGPLOT. 
At the end of arun, the surface impact statistics are summarized and the deposition 


fraction (the mass fraction of released material that has reached the surface in the desired 
simulation time) is computed. 


AGDISP is run from and produces data into files assigned to specific unit numbers. 
These files are the following: 
UNIT DESCRIPTION 

AGDISP input file containing the appropriately constructed data entry cards. 


4 

5 Terminal input. 
6 Screen output. 
7 


CASEFILE default file set of AGDISP input files (accessed with card 0005). 


10 
12 
15 
20 


25 


30 


Binary plot files constructed by AGDISP and plotted by AGPLOT. 
AGDISP.LIS printer listing file. 

WAKE binary plot input file AGWAKE.BIN. 

Binary log file constructed by AGDISP to control plotting by AGPLOT. 
AGPLOT.PLT file constructed by AGPLOT to hold TR instructions. 


AGDISP.ERR and AGPLOT.ERR files generated by AGDISP and AGPLOT if an 
error occurs in the personal computer environment (subsequently checked by 
AGCTRL to detect the error). 


AGPLOT.STP file generated by AGPLOT if a plot is to be displayed and control 
returned to AGPLOT in the personal computer environment. 


Temporary file generated by AGDISP and AGPLOT for intermediate operations. 


The current limits of the AGDISP code are the following: 


1. All of the idealized background mean velocity and turbulence fields are for 
neutral environments. The inputted crosswind velocity profile will be 
interpreted as locally neutral. A Richardson number correction (card 0027) may 
be used for stable/unstable effects of turbulence on diffusion of the material 
from the nozzles. 


2. The total number of nozzles that can be included in the simulation cannot exceed 
60. The error message "Incorrect number of nozzles" will be displayed and 
AGDISP will terminate. 


3. The Betz roll-up cannot enter more than 100 discrete circulation values along the 
wing. The error message "Error in Betz circulation data input” will be 
displayed and AGDISP will terminate. 


4. The Betz roll-up procedure cannot handle more than four discrete vortices on 
each wing. In this case the error message "Betz will roll up more than four 
vortices" will be displayed and AGDISP will terminate. 


5. The plant canopy input data cannot permit more than 20 discrete entries to 
define the canopy vertical profile shape. The error message "Error in plant area 
input" will be displayed and AGDISP will terminate. 


6. The discrete crosswind velocity data cannot permit more than 20 entries to 
define the velocity profile. The error message "Error in crosswind velocity 
input" will be displayed and AGDISP will terminate. 


7. The wide body data cannot permit more than 20 entries to define the fuselage 
cross-sectional area as a function of distance from the aircraft nose. The error 
message "Error in wide body area input" will be displayed and AGDISP will 
terminate. 
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8. The input deck to AGDISP cannot exceed 200 lines. The message "AGDISP 
limited to 200 input cards" will be displayed and AGDISP will terminate. 


‘9. No more than 16 drop sizes may be entered into AGDISP for any run. The 
message "Error in number of drop sizes" will be displayed and AGDISP will 
terminate. 


The list of error messages generated in AGDISP are the following: 


AGDISP does not support card ... 
The input card in question is not a valid entry to AGDISP. 


AGDISP limited to 200 input cards 
A maximum of 200 input cards are permitted in the AGDISP input deck. 


Betz will roll up more than four vortices 
A maximum of four vortices are permitted on a wing semispan. 


Card order inconsistent at card ... 
An error in input has lead AGDISP to expect a card that is not in the input 
deck. 


Error in Betz circulation data input 
The Betz roll-up procedure must have between three and 100 discrete 
locations of circulation data. 


Error in crosswind velocity input 
The user-supplied crosswind velocity profile must have between three and 
20 discrete locations of velocity data. 


Error in number of drop sizes 
The number of drop sizes must be between one and 16. 


Error in number of evaporation equations 
If evaporation is parameterized, there must be as many 0066 cards as 0064 
cards defining the number of drop sizes. 


Error in plant area input 
The plant canopy must have between three and 20 discrete locations of plant 
area fraction data. 


Error in wide body area input 
The wide body cross-sectional area must have between three and 20 discrete 
locations of area data. 


Incorrect number of nozzles 
No more than 60 nozzles may be present in a simulation. 


Initial condition mismatch ; 
Initial conditions were entered for a nonexistent nozzle. 
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Input card numbers not increasing 
All input card identification numbers must increase in number (with the 
exception of card 0000). 


Input does not fully initialize AGDISP 
Input data cards are missing. 


Insufficient data before card ... 
AGDISP is not fully initialized; perhaps a card is missing from the input 
deck or does not have enough information on it. 


Invalid CASE number 
AGDISP could not locate the desired CASE number on card 0005 in 
CASEFILE.INP. 


Premature end of WAKE plot file reached 
The WAKE plot file AGWAKE.BIN ended without supplying enough 
profile information data to AGDISP. 


WAKE plot file extrapolation 
Extrapolation of data in the WAKE plot file AGWAKE.BIN first occurs; 
this is a warning message. 


WAKE plot file (Y,Z) mesh size 
The WAKE plot file AGWAKE.BIN must have between two and 16 mesh 
points in the y and z directions. 


AGPLOT interactively plots AGDISP results on appropriate terminal screens or 
plotting devices. After reading the log file AGPLOT offers the current menu of plotting 
options available to the user (Section 5), who then determines the options to invoke. The 
list of error messages generated in AGPLOT are the following: 


Ground encounter by equivalent gaussian 
The computed equivalent Gaussian intersects the surface; this is a warning 
message. 


Maximum scale less than minimum scale 
The user-supplied scale limits are inconsistent, and AGPLOT requests re- 
entry of the scale data. 


Maximum scale size adjusted 
The user-supplied scale limits are adjusted for consistency with the supplied 
scale delta. 


Option produces nothing to plot 
The invoked option does not produce any plotable results. 


Plot option not available for this LOG file 
AGPLOT cannot invoke a requested option because the data needed for it 
was not saved in AGDISP or is not indicated in the LOG file. 


Scale increment less than zero 


The user-supplied scale increment was inconsistent, and AGPLOT requests 
re-entry of the scale data. 


Too many scale divisions 
The user-supplied scale data forced more than ten scale divisions, and 
AGPLOT requests re-entry of the scale data. 


On the Data General the following error messages are also possible: 


GKS_DRAW error 

GKS_OPEN error 

GKS_Q_WS_WIN error 
GKS graphics errors on the Data General have occurred. AGPLOT 
operation should be reviewed and possibly repeated. 
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4. AGDISP Inputs 


This section of the AGDISP User Manual details the input data to the program. All 
data entry is in free format, with number data separated by commas or blank spaces. This 
convenience offers ease of formatting the data but requires that every data card have all of 
its data values present, even if they are zero. Unless noted below as integer values, data is 
entered as real numbers (with decimal points). The MKS system is used throughout this 
section of the report, even though card 0001 permits English units entry. 


All data cards begin with a four-digit identification in columns 1-4, with the rest of 
the data in free format in column 6 and following. The order of the cards is important and 
must follow ever-increasing identification numbers. AGDISP has been programmed to 
verify this order. In addition, certain available options are inconsistent with each other. 
AGDISP has been programmed to trap these inconsistencies, with the message "Card order 
inconsistent at card ..." followed by the offending card when an error is detected. 
Unsupported cards are flagged with the message "AGDISP does not support card ...", 
while missing data cards are flagged with "Input does not fully initialize AGDISP". With 
free format there exists the chance that all necessary data do not appear on the appropriate 
card; in this case, AGDISP will run out of data cards before it rationalizes all of its data 
pointers. On the other hand, because all of the data is echoed on the terminal screen, error 
messages are traceable in a systematic manner. A careful check should be made the first 
time a new case is started, and at least until the user is confident of program input 
requirements. Only by looking at what the program thinks is inputted will the user be able 
to verify that what was inputted was correct. In almost all cases, AGDISP makes no check 
of input validity, either signs or magnitudes. Every AGDISP run requires the entry, at 
least, of cards 0010, 0020, 0050, 0060 and one or more 0064. 


Detailed description of the currently available data cards follows, including the 
special messages they invoke. 


Caution: All input cards are described below, but for a specific AGDISP run, only 
those cards that are needed to describe the run should be included. Default entry (cards 
0010, 0020, 0050, 0060 and 0064) will be sufficient to describe most runs. For example, 
if the run is made without evaporation, there is no reason to include card 0065. 


In the descriptions that follow, Types are designated CHAR for character data, INT 
for integer data (no decimal point), or REAL for data with decimal point. 


0000 CMNT 
Variable Type Description (Comment Card) 
CMNT CHAR Comment card, any number of which may be placed 


anywhere in the input deck. 


0001 IENG 


Variable Type Description (English Units Card) 
IENG INT Conversion to English units. A default value of 0 (or no 


card) will retain the MKS units. If IENG is nonzero, 
AGDISP will expect weight and thrust in pounds; lengths 
and distances in feet; velocities in miles per hour; and 
temperature in degree F. 


0005 ICASE 


Variable Type Description (Casefile Card) 


ICASE INT Case to run (or modify) found in the default case input 
datafile CASEFILE.INP. All of the input cards found on 
the case file for this case value will be added to the input 
runstream to build a complete input deck for AGDISP. 
Cards present in the input file will replace any cards of 
identical number found in the case file. 


The case file may contain any number of input files to 
AGDISP, each with its input information consistent with 
the discussion in this section. The only attribute that 
separates the case file decks is a first card to each deck 
containing the word CASE in columns 1-4 and the case 
number in column 6 and beyond. With card 0005 present, 
AGDISP will search the default case input datafile to find a 
match of case number; if that match is made, all of the 
subsequent cards up to the next CASE card will be copied 
to the input runstream. If a match is not made, AGDISP 
will print "Invalid CASE number" and terminate. 


0010 TMAX 
Variable Type Description (Program Card) 


TMAX REAL Maximum time in seconds for the simulation to run, a time 
which will be shortened if the material impacts the surface 
or is discarded at the evaporation cutoff (with the proper 
input on card 0065) before TMAX is reached. 


0011 NPTJ NPVX NPVL NPXX 


Variable Type Description (Printer/Plotter Control Card) 
NPTJ INT Modified default printing interval. Under normal operation 


AGDISP output is generated at appropriate intervals to the 
printer and the plot file. Entry of a number greater than 1 
increases the printing interval (the plotting interval cannot 


NPVX 


NPXX 


0015 TA 


Variable 
TA 


Type 
REAL 


be altered). Default value is 1; to turn off printing requires 
an entry of 0. 


Invokes plot file saving of the time histories of the driving 
trajectories, and controls the interval of printing these time 
histories. The driving trajectories are those associated with 
the centers of vortices, the center of the propeller and the 
center of the helicopter rotor disk (whenever any are 
present in a run). Any value for this entry invokes 
additional plot file output, to permit subsequent plotting of 
the driving trajectories by AGPLOT. A number greater 
than 1 increases the printing interval (if nonzero, NPVX 
should equal NPTJ for consistent printed output). Default 
value is 0. 


Controls whether the vertical velocities are saved on the 
plot file for subsequent plotting by AGPLOT. Default 
value is 0 (no saving of vertical velocities); any nonzero 
value entered on this card will force additional plot file 
output. Vertical velocities must be saved if the equivalent 
Gaussian distribution is desired. 


Controls whether the mean velocities are saved on the plot 
file for subsequent plotting by AGPLOT. Default value is 
0 (no saving of mean velocities); any nonzero value 
entered on this card will force additional plot file output. 
Mean velocities must be saved if deposition on objects is 
desired. 


In all cases, additional writes to the plot file increases its 
size. 


Description (Terrain Slope Angle Card) 


Local ground slope in degrees. A positive terrain angle 
raises the nght side of the surface as viewed from behind 
the aircraft, with the origin remaining along the initial 
aircraft vertical centerline. The terrain is characterized by a 
locally straight surface so that all of the simplified flow 
field options in AGDISP remain available. Canopy and 
crosswind effects remain parallel to the tilted surface. 


0020 LMVEL S DIST UO 


Variabl 


LMVEL 


Type 
INT 


Description (Aircraft Characteristics Card) 
Mean velocity flag taking one of seven values: 


5 denotes ground sprayer; 
4 denotes helicopter entry; 
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DIST 


UO 


0021 DZBP 


Variable 


DZBP 


PSBP 


PGBP 


0023. WT 
Variable 


REAL 


PSBP 


Type 
REAL 


Type 


denotes an elliptically loaded, rolled-up tip vortex; 
denotes a rectangularly loaded, rolled-up tip vortex; 
denotes a triangularly loaded, rolled-up tip vortex; 
denotes Betz roll-up from a given circulation, and 
denotes WAKE plot file entry (explained under card 
0050). 


=O — WN WD 


An entry of 4 requires a 0030 card; entry of 3, 2 or 1 
requires a 0023 card; and an entry of 0 requires multiple 
entries of 0025 cards. 


Semispan of the aircraft in meters or feet (also the rotor 
radius for a helicopter). 


Nominal height of the aircraft wing (or nominal height of 
the helicopter rotor blades) about the surface in meters or 
feet (this distance is the assumed nominal nozzle release 
height). For rolled-up tip vortices this height is the initial z 
coordinate of the vortex centerline. For Betz roll-up data 
this height is the z coordinate of the initial vortex sheet . 


Flight speed of the aircraft in m/sec or miles/hour. 


Description (Biplane Characteristics Card) 


Vertical distance in meters or feet from the main wing 
location specified on card 0020 to the second, biplane 
wing. 


Semispan of the biplane wing entered as a fraction of the 
semispan of the main wing (if the wings are equal in length 
this entry would be 1.0). 


Weight carried by the biplane wing entered as a fraction of 
the weight carried by the main wing (for equal weights this 
entry would by 1.0). 


Description (Aircraft Weight Card) 


Weight of the fixed-wing aircraft in Newtons or pounds, 
for a rolled-up tip vortex pair. AGDISP internally corrects 
the proportion of aircraft weight assigned to each wing in a 
biplane configuration (using the card 0021 data). 
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0027 RI 


Variable 


RI 


Type 
REAL 


0028 U Z ZO CA 


Variable 


U 


fe 


Type 
REAL 


Description (Betz Wing Load Distribution Cards) 


Position in meters or feet measured from the wing root 
toward the wing tip. 


Circulation value in m_2/sec or ft2/sec at this location. 


After all 0025 cards are read, the Betz roll-up procedure is 
initialized. The ensuing roll-up invoke additional printer 
output summarizing the roll-up process. As the Betz roll- 
up continues, the material being tracked will be influenced 
by ever-increasing strength vortices, whose positions and 
strengths approach those of rolled-up vortices. AGDISP 
includes the effect of the unrolled-up sheets on the motion 
of the released material. 


Description (Richardson Number Card) 


Richardson number defining the stability of the 
atmosphere: 

cael) 
boa © dz 

au? 


Oz 








For a neutral atmosphere, Ri = 0.0 (default); for stable 
layers, Ri may reach a value of 0.5; for unstable layers, Ri 
may be as or more negative than -1.0. The value of Ri 
modifies the background turbulence level (Bilanin, Teske 
and Hirsh 1978). 


Description (Neutral Crosswind Card) 


Mean wind velocity in the (x,y) plane in m/sec or 
miles/hour. 


Altitude of the mean wind velocity in meters or feet (if a 
canopy is present, this altitude must be above the canopy). 


Surface roughness height z, in meters or feet. 


0029 Z U VV 
Variable Type 
ZL REAL 
U REAL 
V REAL 


0030 WT BDOT 


Variable Type 
WT REAL 
BDOT REAL 


Direction angle in degrees from which the wind is 
blowing; i.e., for 0 deg the wind is a head wind; for 180 
deg the wind is a tail wind; for 90 deg the wind is a 
crosswind from right to left; and for -90 deg (or 270 deg) 
the wind is a crosswind from left to right. 


Description (Discrete Crosswind Cards) 


Position in meters or feet beginning at the surface and 
increasing to the top of the computational domain. 


Axial velocity value in m/sec or miles/hour, positive 
pointing downstream (head wind). 


Crosswind velocity value in m/sec or miles/hour, positive 
pointing in the positive y direction. The square root of the 
sum of the squares of U and V recovers the total velocity 
magnitude at height z. 


After all 0029 cards are read, a calculation is made for the 
turbulence level through the velocity profile. If a canopy is 
present, it is assumed that the entries on cards 0029 have 
been modified by the user to reflect the influence of the 
canopy on the crosswind velocity profile and its 
subsequent turbulence calculation. 


Description (Helicopter Card) 

Weight of the helicopter in Newtons or pounds. 

Blade rotation rate in rpm. 

The helicopter is idealized as a dynamic transition from a 


rotor downwash field to a rectangularly loaded, rolled-up 
vortex pair. 


0035 THT RPRP DZ XPRP 


Variable Type 
fase REAL 
RPRP REAL 
DZ REAL 


Description (Jet Engine Card) 

Thrust of the engine in Newtons or pounds. 

Exit radius of the engine in meters or feet. 

Incremental distance in meters or feet of the engine 


centerline above or below the nominal release height given 
on card 0020 (the jet engine is assumed to be at the aircraft 
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XPRP REAL 


centerline, y = 0). The placement of the jet engine is 
superseded by the presence of card 0045. 


Axial distance from the trailing edge of the wing to the exit 
plane of the jet engine in meters or feet, where positive is 
measured upstream. 


0040 CD AS ETA TDOT RPRP DZ XPRP 


Variable Type 

CD REAL 
AS REAL 
ETA REAL 
TDOT REAL 
RPRP REAL 
DZ REAL 
XPRP REAL 


Description (Propeller Card) 
Drag coefficient of the aircraft (typically 0.1). 


2 2 


Airplane planform area in m~ or ft“. 


Propeller efficiency (typically 0.8 to 0.9). 
Shaft rpm. 
Propeller blade radius in meters or feet. 


Incremental distance in meters or feet of the shaft centerline 
above or below the nominal release height given on card 
0020 (the propeller is assumed to be at the airplane 
centerline, y = 0). The placement of the propeller is 
superseded by the presence of card 0045. 


Axial distance from the trailing edge of the wing to the 
propeller blade plane in meters or feet, where positive is 
measured upstream. 


0045 NPRP (TV(N), N = 1 TO NPRP) 


Variable Type 
NPRP INT 
TV REAL 


Description (Powerplant Placement Card) 


Number of engines (1 to 4, although 1 is the default on 
cards 0035 and 0040). 


An equal number of dimensions following the pattern 


for one engine: 
incremental distance (meters or feet) at the aircraft 
centerline above or below the nominal height; 


for two engines: 
positive horizontal position (meters or feet) of the 
engine from the aircraft centerline; and 
incremental distance (meters or feet) of the engine 
above or below the nominal height (the other engine is 
symmetrically positioned); 
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for three engines: 
incremental distance (meters or feet) at the aircraft 
centerline above or below the nominal height for the 
engine at the aircraft centerline; 
positive horizontal position (meters or feet) of the 
second engine from the aircraft centerline; and 
incremental distance (meters or feet) of the engine 
above or below the nominal height (the third engine is 
symmetrically positioned); 


for four engines: 
positive horizontal position (meters or feet) of the first 
engine (closer to the aircraft centerline) from the 
aircraft centerline; 
incremental distance (meters or feet) of the engine 
above or below the nominal height (the second engine 
1s symmetrically positioned); 
positive horizontal position (meters or feet) of the third 
engine (farther from the aircraft centerline) from the 
aircraft centerline; and 
incremental distance (meters or feet) of the engine 
above or below the nominal height (the fourth engine is 
symmetrically positioned). 


0047 (IPRPS(N), N = 1 TO NPRP) 


Variable Type Description (Propeller Rotation Card) 
IPRPS INT Up to four integers (depending on the number of engines 


in the simulation). Propellers are numbered from the left 
wingtip to the right, looking upstream; the integer is 1 if 
the propeller rotates clockwise, and -1 if the propeller 
rotates counterclockwise. Defaults are as follows: 


1 propeller 1 
2propellers 1 -1l 
3 propellers 1 1 -1 
4propellers 1 ere bea 
0050 LQQSE QQMxX 
Variable Type Description (Turbulence Card) 
LQQSE INT Turbulence flag carrying one of the following values: 


1 invokes superequilibrium turbulence; 
0 assumes a fixed value of turbulence given on this 
card; and 
-1 specifies the attached WAKE plot file (invoked with 
a -1 entry on card 0020). 
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QQMX 


REAL 


Maximum value of the background turbulence q2 in 


m2/sec~ or ft2/sec2 (this value is incremented by the 
presence of a crosswind). The turbulence level generated 
by a crosswind flow field is generally sufficient to 
represent the ambient turbulence level (the input on card 


0050 would then be 0.0 for q?). However, in some cases, 
it may be necessary to augment this value for specific local 
atmospheric conditions. A short guide to the selection of 
appropriate background turbulence levels is offered in 
Section 9. 


With no crosswind, QQMX must be nonzero on this card. 


A -1 on card 0020 invokes the WAKE plot file and must 
be accompanied by a -1 on this card. The WAKE plot file 
(AGWAKE.BIN) contains the crossplane velocities V and 


W, and the turbulence q?. 


The sequential binary WAKE plot file is constructed as 
follows (with FORM='BINARY' on personal computers): 


Record 1: The number of y (or horizontal) mesh points and 
z (or vertical) mesh points in the plot file (two integers). 
The code restricts both of these entries to 16 or less; larger 
values invoke the error "WAKE plot file (Y,Z) mesh size"; 


Record 2: The y mesh values in meters; 
Record 3: The z mesh values in meters; 


Record 4 and following: The profile data follows, with 
each time slot repeating the same pattern. It begins with a 
single record which is the time saved in seconds. The data 
for each variable on the file (specified in the order V, W, 


q2) follows, by giving all of the y values of V at the first z 
position, then all of the y values of V at the second z 
position, on to the last z position; then on to all of the y 
values of W at the first z position, etc., until the values of 
all of the variables at all of the y and z positions have been 
given. The next time value follows, repeated to the end of 
the plot file. Velocities are in m/sec and turbulence in 


m2/sec*. The warning message "WAKE plot file 
extrapolation" is output the first time spatial extrapolation 
must be used by the AGDISP code during interpolation for 
the variables contained on the WAKE plot file. For times 
beyond the entries in the WAKE plot file, the spatial 
profiles nearest in time to the time being solved in 
AGDISP will be used. 


Before integrating the equations, AGDISP reads the entire 


WAKE plot file through to its end. If the file is 
constructed incorrectly, the error "Premature end of 
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0055 Z AA 


Variable 
Z 


AA 


Weoge CLEP 


Variable 


EEF 


0060 NVAR 
Variable 


NVAR 


DZ 


XO 


DENF 
FLOW 


Type 
REAL 


Type 


REAL 


WAKE plot file reached" is invoked, and AGDISP 
terminates. 


Description (Canopy Cards) 


Position in meters or feet beginning at the surface and 
increasing to the top of the canopy. 


Plant area fraction corresponding to the z position. 


After all 0055 cards are read, the canopy calculation is 
initialized, computing the displacement thickness of the 
canopy and forcing modification to the crosswind velocity 
and turbulence within the canopy. During the AGDISP 
run, trajectories of the aircraft vortices will be altered upon 
entering the canopy, and the position of the helicopter 
dividing streamline will also be changed. 


Description (Canopy Collection Efficiency Card) 


Canopy capture efficiency B (Eq. 47). Default value is 1. 


DZ XO DENF FLOW 


Type 
INT 


Description (Nozzle Card) 


Number of nozzles in the entire calculation (this number 
may be odd or even). Unless followed by cards 0061, 
AGDISP will position the nozzles uniformly along the 
wing (or rotor centerline) between minus and plus three- 
fourths of the semispan (or rotor radius), numbered from 
left to right from behind the aircraft. 


Vertical position in meters or feet off-setting the nozzles 
from the height of the wing (or rotor plane) given on card 
0020. 

Nominal axial position of the nozzles relative to the trailing 
edge of the wing (or the shaft centerline of the helicopter) 
in meters or feet, positive measured upstream. 

Specific gravity of the released material. 


Flow rate from the spray system in gallons/minute. 


0061 Y Z 


Variable 
Y, 
Z: 


Type 
REAL 
REAL 


Description (Discrete Nozzle Location Cards) 
Y value in meters or feet along the spray boom. 


Vertical position in meters or feet off-setting the release 
point from the height of the wing (or rotor plane) given on 
card 0020. The vertical position offset on card 0060 is 
disregarded. Discrete nozzle locations permit fine tuning 
of nozzle locations, and are necessary if nozzles are not 
uniformly distributed along the wing or if the spray boom 
is not parallel to the surface. The number of 0061 cards 
must equal the number of nozzles NVAR on card 0060. 


OD62°N.aU aVaeiW eXSVS 


Variable 
N 
U 


XS 


VS 


0064 DIAM 
Variable 
DIAM 

FRAC 


Type 
INT 


FRAC 


Type 


Description (Nozzle Initial Condition Cards) 
Nozzle number, recorded from the left nozzle position. 


Initial axial velocity in m/sec or miles/hour at nozzle 
number N, positive measured downstream. 


Initial horizontal velocity in m/sec or miles/hour at nozzle 
N. 


Initial vertical velocity in m/sec or miles/hour at nozzle N. 


Initial spatial variance of the material path in m2 or ft? at 
nozzle N. 


Initial velocity variance of the material in m2/sec2 


ft2/sec2 at nozzle N. 


Or 


All initial conditions are set to zero unless modified by card 
0062 for each nozzle. Card 0062 is ideal for an AGDISP 
run with few nozzles, such as a ground sprayer 
simulation. 


Description (Drop Size Card) 
Initial diameter of the drop in microns. 
Mass fraction assigned to this drop size. The sum of the 


mass fractions for all drop sizes in the simulation must 
equal 1.0 or less. 


0065 DTEMP VFRAC 


Variable Type 
DTEMP REAL 
VFRAC REAL 


0066 DEA DEB DEC 


Variable Type 

DEA REAL 
DEB REAL 
DEC REAL 


Up to 16 drop sizes may be specified in one AGDISP run. 
Drop sizes and their mass fractions may be estimated for 
common spray materials using DROPSIZE (Teske 1990b) 
and combining categories with SDC (Teske 1989d). 


. . 


Description (Evaporation Card) 


Wet bulb temperature difference in deg C or deg F. The 
program WETBULB (described in Section 8) may be used 
to generate this input. 


Volatile fraction of the released material (the fraction of 
material that will evaporate; if no material evaporates, then 
VFRAC = 0.0). 


If the wet bulb temperature difference entered on this card 
is zero, AGDISP assumes that evaporation will be 
parameterized, and expects subsequent cards 0066. If 
VFRAC is negative, material is removed from the 
simulation when it evaporates down to the absolute value 
of VFRAC. 


Description (Evaporation Parameterization Card) 


Coefficient A in microns in the material diameter equation 


Diameter = A + Bt + Ct? (50) 
Coefficient B in microns/sec. 


Coefficient C in microns/sec2. 


These coefficients are compatible with the evaporation 
model in FSCBG (Teske and Curbishley 1991), and 
produce consistent evaporation results between AGDISP 
and FSCBG. The coefficient A must be identical to each 
material diameter entered on cards 0064 and in the same 
order (the order is NOT checked by AGDISP). The 
number of 0066 cards must equal the number of drop sizes 
(0064 cards). 


0070 Z 


Variable 
Th 


0075 GDK 
Variable 


GDK 


0080 XBOD 
Variable 


XBOD 


DZ 


0081 XxX SS 


Variable 
Xx 


Type 
REAL 


Type 


DZ 


Type 


Type 


Description (Apparent Surface Height Card) 


Apparent height above the surface (in meters or feet) at 
which the material will deposit. The default condition 
requires impact of material at the surface (with an apparent 
height equal to zero). For a nonzero entry, when material 
intersects the apparent height, it is removed from the 
computation, just as though it had impacted the surface. A 
subsequent ground deposition plot with AGPLOT will 
recover the deposition of the material at the apparent height 
entered on this card. 


Description (Circulation Decay Card) 


Factor (nondimensional) that forces a time reduction in the 
circulation strength of the vortices. Default value is zero 
(no decay). If this entry is negative, its absolute value is 
taken as the circulation decay constant multiplied by the 
ambient turbulence level (m/sec or ft/sec). Teske 1988b 
suggests that for small material sizes, an appropriate entry 
on card 0075 is -0.56 m/sec. 


Description (Wide Body Setup Card) 


Position in meters or feet of the front of the fuselage from 
the trailing edge of the wing (or the shaft centerline of the 
helicopter), positive measured upstream (by convention, 
then, this number must be positive). 


Vertical distance in meters or feet off-setting the nose 
centerline of the fuselage from the nominal height of the 
wing (or rotor plane) given on card 0020. The fuselage is 
assumed to be axisymmetric. If card 0080 is present, 
cards 0081 must also be present to define the fuselage 
body shape. 


Description (Wide Body Effect Cards) 


Position in meters or feet, measured positive from the nose 
of the fuselage to the tail, increasing toward the tail. 


SS 


0090 N2FRO 
Variable 


N2FRO 


Type 
INT 


Cross-sectional area of the fuselage in m2 or ft2, 
corresponding to the x position. 


The first 0081 card should contain 0.0 for x and 0.0 for 
the area. The last 0081 card should contain 0.0 for the 
area. After the final 0081 card is read, wide-body effects 
are initialized by a source-sink velocity profile near the 
fuselage and an axisymmetric turbulent wake decaying 
behind the fuselage. 


Description (To And Fro Card) 


AGDISP always computes results for the "to" direction, 
with the aircraft flying away from the observer; any 
nonzero entry on this card will also invoke a 
corresponding "fro" direction, with the aircraft flying 
toward the observer. 


5. AGPLOT Inputs 


The AGPLOT program may be invoked to plot the resulting trajectories and 
deposition pattern from an AGDISP run. The program asks questions interactively of the 
user, and with this information processes the data and plots the results. 


At the beginning of AGPLOT, the log file is read, the input data deck used to 
generate the file is written to the screen (to serve as verification of the proper run to plot), 
and the log file is scanned for maximum values and available options. With this 
information in hand, AGPLOT then displays its current menu driver: 


0 EXIT AGPLOT 
1 Mean drop trajectories 
2 Mean + standard deviation trajectories 
3. __-Vortices/helicopter/engine centroids 
4 Continuous ground deposition 
5 Volume median diameter 
6 Number median diameter 
7 Layered canopy deposition 
8 Total canopy deposition 
9 Vertical flux profile 
10 Equivalent gaussian distribution 
11 Crosswind velocity profile 
12 Canopy plant area fraction profile 
13. Drop diameter time history 
14 Drop canopy parameter time history 
15 Drop velocity time history 
16 Object deposition time history 
17 _—_‘ Drift fraction time history 
18 Coefficient of variation profile 
19 Swath deposition overlap pattern 


followed by the request "Enter option to run". Only those options available from the saved 
datafiles will be displayed. Each of the options is described as follows: 


0 is entered to exit from AGPLOT. Since several menu options may be 
invoked for one AGDISP run, 0 is the entry used to complete the plotting 
for a particular log file. 


1 is entered to plot the mean trajectories. If more than one drop size was 
solved (with more than one card 0064), AGPLOT will display the available 
BIN file numbers ("LOG file contains more than one BIN file"), and their 
corresponding initial drop diameters, and asks 


Enter BIN file to access 


to select the datafile to read. 


is entered to plot the mean trajectories (as solid curves) and the standard 
deviation of the mean trajectories (as dashed curves). 


is entered to plot the trajectories of the vortex centroids (as solid curves) and 
the center(s) of the propeller, jet engine or helicopter (as dashed curves), if 
applicable. This option is available only with an appropriate nonzero entry 
on card 0011 in AGDISP. 


is entered to compute and plot the composite ground deposition pattern with 
continuous deposition following Eq. (46). Up to 16 drop sizes may be 
combined in this plot, as specified by cards 0064 in the input deck. If "to" 
and "fro" calculations have been invoked with card 0090, the user is asked 


LOG file contains To and Fro calculations ... 
Which direction to access (T/F) 


with the "to" direction (T) as default. Seven scales are available in 
AGPLOT: 


normalized 
liters/hectare 
fluid ounces/acre 
gallons/acre 
mg/sq m 

nl/sq cm 
drops/sq cm 


ANNA RWNH 


Also, if the volatile fraction is nonzero (meaning there is evaporation), 
AGPLOT asks 


Plot nonvolatile contribution only (N/Y) 


with no (N) as the default. If the nonvolatile contribution is invoked, 
"NONVOLATILE" appears in the upper right corner of the plot. For the 
normalized deposition plot, an integral under the deposition curve reflects 
the amount of material that has impacted the surface. This number appears 
in the lower left corner of the plot. 


When the calculation is completed, AGPLOT asks 
Save ASCII data (N/Y) 


to permit the user to port the results into an alternate graphics package. If 
the response is yes (Y), AGPLOT asks | 


Enter FILENAME.EXT 


for the user to enter an appropriate filename designation to save the results. 
Results are saved in a two-column format file, where the first column is the 
horizontal position relative to the aircraft centerline (the y distance), in 
meters measured from left to right; and the second column is the 
corresponding values in the scale previously selected. 
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10 


is entered to compute and plot the total volume median diameter of all 
material deposited on the ground. The average VMD over the deposition 
appears in the lower left corner of the plot. 


is entered to compute and plot the total number median diameter of all 
material deposited on the ground. The average NMD over the deposition 
appears in the lower left corner of the plot. 


is entered to compute and plot contours of the composite canopy deposition 
pattern following incremental Gaussian deposition (Eq. 44) through the 
canopy. Deposition is computed in layers through the canopy. The user is 
then quizzed to give specific contour values to be plotted, or to permit 
AGPLOT to autoscale the contour values. Option 7 must be exercised to 
generate the profile distribution needed for option 8. The maximum 
deposition value appear in the lower left corner of the plot. 


is entered to plot the total canopy deposition determined in option 7. 


is entered to compute and plot the total flux through a vertical plane defined 
by the user in response to the question 


Enter Y position of plane (m) 


Total vertical flux is the sum over time of the material passing through the 
specified plane. The y position is shown in the upper left corner of the plot. 


is entered to compute and display the equivalent Gaussian distribution. 
Since the position and standard deviation of all material in the simulation is 
computed by AGDISP as a function of time after release, the equivalent 
mean position and standard deviation may be determined by appropriate 
summation and integration at each time over all material in the simulation. A 
measure of the compatibility of the equivalent Gaussian distribution is made 
by computing a figure of merit ranging from 0 to 1. When the equivalent 
Gaussian no where represents the multiple-drop distribution, the figure of 
merit equals zero. When the equivalent Gaussian is everywhere identical to 
the multiple distribution, the figure of merit is unity. Equivalent Gaussian 
calculations are made for every time in the designated plot file, with all times 
given a sequence number. The critical sequence number, the time 
associated with each of them, and their respective figure of merit and 
terminal velocity criterion are displayed on the screen whenever: the first 
sequence number is reached; the last sequence number is reached; the figure 
of merit reaches a local minimum; the figure of merit reaches a local 


maximum; or material first comes within a standard deviation of the ground. 


The terminal velocity for the material size examined is computed assuming 
no evaporative effects and a specific gravity of unity. It is displayed with 
the message "Terminal velocity" as reference to the user. The terminal 
velocity criterion is computed as the average over all material of their vertical 
velocities relative to their terminal velocity (if all material were moving at 


aie) 


11 
12 


13 


14 


15 


their terminal velocity, the WCRIT value would be 1.0). Option 10 requires 
that vertical velocities be saved on the plot file by the appropriate input on 
card 0011 in AGDISP. Material whose vertical velocities are not within 
fifty percent of their terminal velocity will not contribute to the figure of 
merit calculations at that time step. 


After processing the entire plot file and displaying the above summary 
information, the plot file is reread to locate the position of maximum figure 
of merit. At this point the equivalent Gaussian data is displayed on the 
screen. If the equivalent Gaussian vertical standard deviation contacts the 
surface, the warning message "Ground encounter by equivalent gaussian" 
will be displayed on the screen. 


The resulting plot represents each released material distribution as dashed 
curves and the single equivalent Gaussian distribution as solid curves. 
Around each center point of the material contributing to the equivalent 
Gaussian, normalized contour lines are plotted at one-fourth and one times 
the standard deviation. Material not contributing to the equivalent Gaussian 
are plotted with much smaller ellipses to aid in identifying them. The figure 
of merit appears in the lower left corner of the plot. 


is entered to plot the crosswind velocity profile (if available with a 0028 
card or 0029 cards in AGDISP). 


is entered to plot the canopy plant area fraction profile (if available with 
0055 cards in AGDISP). 


is entered to plot the time history of a drop diameter (if evaporation iS 
occurring). If there are more than one nozzle in the simulation, AGPLOT 
requests 


Enter nozzle to process 
The selected nozzle number is displayed in the lower left corner of the plot. 


is entered to plot the time history of the canopy deposition of a drop (if 
canopy cards 0055 are present in AGDISP). 


is entered to plot the velocity time history of a drop (if the mean velocity 
data is stored on the plot file with card 0011 in AGDISP). AGPLOT asks 


1 U axial velocity 
2  V horizontal velocity 
3. W vertical velocity 


to enable the user to select one of the velocities for the plot. The selected 


velocity (U AXIAL, V HORIZONTAL or W VERTICAL) is written to the 
lower right corner of the plot. 
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16 


17 


18 


is entered to compute and plot the time history of the deposition on a 
specified collector geometry at a specified position in the AGDISP solution 
field (mean velocity data must be stored on the plot file with card 0011 in 
AGDISP). The collector must first be placed in the solution field by 
responding to the request 


Enter (Y,Z) collector location (m) 


with the y and z values entered in meters. Next, AGPLOT asks for the type 
of collector; currently, the options are: 


1 Card 
2 Cylinder 
3 Sphere 
Then AGPLOT requests 


Enter significant target dimension (m) 


to give a typical length of the collector in meters. Finally, AGPLOT needs 
to know the orientation of the collector relative to the (x,y,z) coordinate 
system, and asks 


Enter (X,Y,Z) normal vector 


Here all three components must be entered; AGDISP renormalizes into a 
unit vector. A preferred normal direction must be supplied, even for 
collectors that do not appear to have planar geometry. If the collector were 
placed parallel to the surface, the normal vector would be (0,0,1). If the 
collector were vertical and deposition on the left side were needed, the 
normal vector would be (0,-1,0). 


The total amount of material collected on the collector appears in the lower 
left corner of the plot. The collection efficiency appears in the lower right 
corner of the plot. 


is entered to compute and plot the time history of the drift of all of the drops 
in the simulation. The final drift fraction appears in the lower left corner of 
the completed plot. 


is entered to compute and plot the coefficient of variation as a function of 
lane separation between flight lines of the aircraft. Option 4 (continuous 
ground deposition) must be invoked to generate the ground deposition 
before option 18 can be invoked to overlap the deposition to compute lane 
separation effects. Number depositions (in drops/sq cm) should not be 
used for lane separation computations. If both "to" and "fro" ground 
depositions are available, the user will be asked to select one of the options 


Which direction to access (T/F/C) 


where the default is "to" (T), a response of F will invoke "fro" and a 
response of C will invoke a combined deposition, alternating between the 


2a! 


two saved depositions to compute a composite pattern. The complete plot 
will include, in the lower left corner, the value of lane separation at a 
coefficient of variation of 0.3, or at the curve minimum location. This value 
of COV is suggested in the literature (Teske, Twardus and Ekblad 1990). If 
the resultant curve for COV is not "smooth" enough, a smaller horizontal 
distance must be chosen in option 4 for continuous ground deposition. 


19 is entered to compute and plot the effective swath pattern if multiple flight 
lines are overlapped at a lane separation specified by the user 


Enter lane separation distance (m) 


The lane separation (from option 18), or any other separation value, may be 
entered here. Option 4 (continuous ground deposition) must be invoked to 
generate the ground deposition before option 19 can be invoked to overlap 
the deposition. If both "to" and "fro" ground depositions are available, the 
user will be asked to select one of the options 


Which direction to access (T/F/C) 


The complete plot will include, in the lower left corner, the average 
deposition value over the swath width (in the user specified deposition 
units). 


Where appropriate, the following requests are made by AGPLOT: 

Autoscale canopy contours (Y/N) 

AGPLOT contour plots appropriate deposition levels through the canopy in 
option 7. If N is the response, AGPLOT gives the maximum deposition 
value and asks for contour levels (and will plot them) until a negative 
deposition value is entered. 

Autoscale computation grid (Y/N) ; 
AGPLOT computes appropriate solution intervals, if invoked, for options 
that require a specified grid before computations can occur. If N is the 
response, AGPLOT asks 

Enter Y scale min, max 
or 
Enter Z scale min, max (in option 9 only) 

Autoscale plot axes (Y/N) oa 
AGPLOT computes appropriate plotting scales if invoked. If N is the 
response, AGPLOT asks 


Enter Y scale min, max, incr 


and 
Enter Z scale min, max, incr 
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AGPLOT has a scale checking routine and invokes a warning message 
whenever the scale delta is not an integer fraction of the overall scale size. 
Additionally, three errors will be trapped: "Maximum scale less than 
minimum scale"; "Scale increment less than zero"; and "Too many scale 
divisions" whenever more than ten scale divisions are needed. 


Available scales 

normalized 

liters/hectare 

fluid ounces/acre 

gallons/acre 

mg/sq m 

nl/sq cm 

drops/sq cm 

Enter scale to use 
One of the available units must be selected by the user. The "normalized" 
scale produces a nondimensional result. 


ANN BWNH 


Collector types 
1 Card 
2 Cylinder 
3. “Sphere 


Enter type to use 
Option 16 requires the selection of a collector type. 


Enter FILENAME.EXT 
At this prompt the user enters the name of the ASCII datafile into which the 
results are stored. 


Enter lane separation distance (m) 
Option 19 requests the distance between flight lines to construct the 
overlapped swath deposition pattern. 


Enter nozzle to process 
Several options require the selection of one of the nozzles in the simulation. 


Enter plot title 
The user enters the desired title for the current plot (a maximum of 32 
characters). 


Enter significant target dimension (m) 
Option 16 requires the size of the collector in meters. 


Enter (X,Y,Z) normal vector 
Option 16 requires the orientation of the collector relative to the coordinate 
system of the AGDISP solution. 


Enter Y position of plane (m) 
Option 9 requires a y position in the AGDISP solution through which to 
evaluate vertical flux. 


Enter (Y,Z) collector location (m) 


Option 16 requires the placement of a collector (in meters) in the AGDISP 
solution. 


Sal 


LOG file contains more than one BIN file ... 

Enter BIN file to access 
Several options require the selection of one drop size. When more than one 
are solved by AGDISP, the user is asked to select the one to access. 


LOG file contains To and Fro calculations ... 

Which direction to access (T/F) 
Several options require the selection of aircraft direction (if the "fro" or F 
option is available with card 0090 in AGDISP). 


Plot nonvolatile contribution only (N/Y) 
When evaporation is active (with card 0065 in AGDISP), AGPLOT permits 
the user to plot the nonvolatile (nonevaporative) part of the material alone 
(with a response of Y), or to plot the total material. 


Save ASCII data (N/Y) 
Several options permit the results to be written to a readable file for porting 
to other graphics applications. 


Tag location at time increments (N/Y) 
This option is asked for mean trajectories, to tag released material positions 
on the plot as a function of time. The size of the tag is scaled to the size of 
the material standard deviation. If the response is Y, AGPLOT asks 


Enter tag time increment 
to enter the incremental seconds between tags. 


To and Fro directions have been saved ... 

Which direction to access (T/F) 
If the "fro" direction is active (with card 0090 in AGDISP) and both the "to" 
and "fro" canopy depositions (option 7) have been invoked, the user may 
select "to" or "fro" for option 8. 


To and Fro directions have been saved ... 

Which direction to access (T/F/C) 
If the "fro" direction is active (with card 0090 in AGDISP) and both the "to" 
and "fro" ground depositions (option 4) have been invoked, the user may 
select "to" or "fro", or an alternating deposition pattern "combined" for 
option 18 (if both "to" and "fro" have been saved with the same deposition 
units and either total or nonvolatile contributions). 


Velocities 
1 U axial velocity 
2 V horizontal velocity 
3. W vertical velocity 
Enter velocity to plot 
In option 15 one of the three velocities (U,V,W) must be specified. 
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On the Data General the following additional messages are possible: 


Enter MASTER filename | 
The filename containing the data to be read (12 characters or less without 
extension) is entered here. 


Enter METAFILE.EXT for this plot : 
If metafiles are open (from a previous query), then each plot must be 
assigned a name and extension for storing on the DG system in the user 
directory. A metafile extension of GKM is recommended for entry into 
CEO. 


Enter printer device type 
Response to this query determines the hardcopy device used. A nonblank 
entry will invoke the QUE name request. A blank entry will prevent any 
output from being hardcopied. 


Enter QUE name code 
A hardcopy unit must be specified by its queue name within the DG 
environment. If a blank name is given, hardcopy of plots will not occur. 


Enter workstation type 
Response to this query determines the terminal screen used. 


Open output METAFILES (N/Y) 
Plots may be spooled to a file for later manipulation. This option permits 
that operation. 


Plot completed: press RETURN to continue 
The DG system blanks the screen after completing a plot unless a pause is 
programmed into AGPLOT. 


Additionally, on the DG system the pressing of any function key at a query line 
returns the user to the menu driver. 


a 


6. Test Cases 


___ Four test cases are included in the User Manual to illustrate the variety of data entry 
available with AGDISP. These examples are included to demonstrate the structure of the 
input data and the data entry requirements of each type of data card and option. They are 
not meant to suggest that certain options go with certain other options, nor are all options 


covered. 


The four examples illustrate the following features: 


1e 


The first test case examines the wake behind an F-15 fighter flying at Mach 
0.5. Because of the high speed, Betz roll-up is invoked, with a wing load 
distribution obtained from vortex-lattice theory. Figure 6-1 displays the 
appropriate input deck. Figure 6-2 shows the mean trajectories and Figure 6- 
3 shows the ground deposition in drops per square centimeter. 


The second test case examines the wake behind a wide body C-130 transport 
plane. Figure 6-4 displays the appropriate input deck. Figure 6-5 shows the 
mean trajectories; Figure 6-6 includes the spread of released material due to 
turbulence; Figure 6-7 shows the paths of the tip vortices and propellers; and 
Figure 6-8 shows the ground deposition in ounces per acre. 


The third test case examines the wake behind a Cessna 188 Ag Truck in 
forward flight (Mission test run 14, droplet size 171 microns, Teske 1988). 
Figure 6-9 displays the input deck while Figure 6-10 shows the mean 
trajectories and Figures 6-11 and 6-12 show the ground deposition for all of 
the material (in liters per hectare), and for the nonvolatile part (normalized). 
Figures 6-13 and 6-14 show the VMD and NMD plots, while Figure 6-15 
gives a vertical flux plot (in liters per hectare) at a y plane position of -10 
meters in the AGDISP coordinate system. Figure 6-16 illustrates the 
equivalent Gaussian distribution; Figure 6-17 shows the crosswind velocity 
profile; Figure 6-18 shows the time history of the evaporating diameter from 
one nozzle in the simulation; Figures 6-19, 6-20 and 6-21 give the 
corresponding drop velocity time histories; and Figure 6-22 shows the 
deposition on a spherical collector. Figures 6-23 and 6-24 illustrate the drift 
fraction for all of the material, and for the nonvolatile part. Finally, Figures 
6-25 and 6-26 show the coefficient of variation for this example case, and the 
resulting overlapped swath pattern. 


The fourth test case examines a Hiller 12E helicopter (Chico test run B-2, 
Teske 1989c). Figure 6-27 displays the input file, while Figure 6-28 shows 
the mean trajectories and Figure 6-29 shows the ground deposition in gallons 
per acre. Figure 6-30 displays contours of canopy deposition, while Figure 
6-31 summarizes the total deposition through the canopy, on the same scale as 
the ground deposition. Figure 6-32 illustrates the crosswind velocity profile; 
Figure 6-33 illustrates the canopy plant area fraction profile; Figure 6-34 gives 
the canopy deposition fraction from one nozzle in the simulation; and Figure 
6-35 summarizes the drift fraction for the simulation. 
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Figure 6-1 AGDISP input file for Example Case 1. 
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Figure 6-4 AGDISP input file for Example Case 2. 
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Figure 6-9 AGDISP input file for Example Case 3. 
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Figure 6-27 AGDISP input file for Example Case 4. 
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7. AGDISP Operation on Mainframe Computers 


The latest versions of AGDISP and AGPLOT are operational on two mainframe 
computer systems. 


DATA GENERAL MV-15000 
AGDISP and AGPLOT are run on the Data General Forest Service network by the 
commands 
XEQ AGDISP 


XEQ AGPLOT 


with resolution of the filenames and extensions performed within the programs. Both 
programs will ask 


Enter MASTER filename 
to enter the filename of the input file to the programs. 
The plotting package GKS is invoked to plot AGPLOT results onto user-selected 
graphics output devices. 


DIGITAL EQUIPMENT VAX 


AGDISP and AGPLOT may both be invoked by command files of the form shown 
on Table 7-1. Input on a file named FILENAME.INP will produce a binary log file 
FILENAME.LOG that controls all saved plot information. The printer file AGDISP.LIS is 
created during AGDISP operation. 


Plotting is accomplished by calls directly to a Tektronix 4025, saving these calls in 
an ASCIl-readable file AGPLOT.PLT for later display by site-specific graphics programs. 


TABLE 7-1 


Command Procedures for AGDISP on VAX 


To Invoke AGDISP: @AGDISP FILENAME 
will generate 
SASSIGN/NOLOG FILENAME. INP 
SASSIGN/NOLOG CASEFILE.INP 
SASSIGN/NOLOG © AGDISP .LIS 
SRUN AGDISP 
To Invoke AGPLOT: @AGPLOT FILENAME 
will generate 
SASSIGN/NOLOG FILENAME.LOG 
SASSIGN/NOLOG AGPLOT.PLT 


SASSIGN/USER SYSSCOMMAND 
SRUN AGPLOT 
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FORO04 
FOROO7 
FOROO9 


FORO12 
FOROLS 
SYSSINPUT 


8. AGDISP Operation on Personal Computers 


A series of additional programs have been written to make the use of AGDISP and 
AGPLOT more straightforward on personal computers. These programs are described 
below and illustrated in the flowchart shown in Figure 8-1. 


AGEDIT 


AGEDIT is a simple full-screen editor that permits the set up and modification of 
AGDISP input files. It contains a built-in Help facility that briefly explains every input 
card type currently available in AGDISP. 
AGEDIT is invoked with the command 
AGEDIT filename 

If a file name is not entered on the command line, AGEDIT asks 
Enter AGDISP Input File To Edit 

The Help facility must be available to AGEDIT or the message is 
AGHELP.DOC Cannot Be Found Or Opened 


If the INP file given the program cannot be found, AGEDIT assumes the file is to 
be created rather than modified. 


Active function keys are the following: 

Fl Switches to the Help file and back, displaying information about the 
current card pointer. Once in Help, movement is up (back to screen 
1 -- the information screen) or down (to the last screen) with the up 
and down cursor arrow keys. 

F2 Switches to the Menu options screen. 

While editing, a brief Help line is always present on the screen; this line displays 
any identifiable card number and the names/number of variables that should be on the card. 
If the card is not recognizable, AGEDIT displays 

Unsupported AGDISP Card 
The editing options available are the following: 


1 The cursor keys (arrow keys) allow movement up, down, left and right in the 
file. 


2 Typing any realistic characters are added at the current cursor position. 


6 
oi 


INSERT may be toggled on (to insert the characters typed) or off (to overwrite 
with the characters typed). 


DELETE or BACKSPACE deletes the character to the left of the present cursor 
position. 


ENTER generates a new line in the file. 
HOME moves to the first character in the current line. 


END moves to the last character in the current line. 


The menu options available are the following: 


0 


Exit/Edit Without Saving Changes 
AGEDIT always verifies program exit with 
Exit (N/Y) 
The response must be Yes; otherwise AGEDIT returns to editing. 
Save Results In Existing File And Exit 
AGEDIT saves any editing changes in the current input file and exits. 
Save Results In New File And Exit 
AGEDIT asks 
Enter New AGDISP Input File 
and then saves the editing session into this file (with the extension INP) and 
exits. Once back to AGCTRL, the master input file may be changed to the new 
AGDISP input file before running AGDISP. 
Recover Aircraft Characteristics From Database 
AGEDIT contains access logic to FSCBG.ACL, the complete aircraft library 
from Hardy 1987 corrected by B. Thompson and B. Kirk (private 


communication). If the binary data file cannot be located, AGEDIT writes 


FSCBG.ACL Cannot Be Found Or Opened. 
Press Any Key To Continue ... 


after which AGEDIT returns to editing. 


Five aircraft types are shown on the screen at a time. Movement through the 
aircraft library is with the up and down cursor arrow keys. The selection of a 
specific aircraft type is made by positioning the highlight line on the desired 
aircraft type and pressing the ENTER key. The menu option may be aborted 
with the F1 function key. 
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When an aircraft type is selected, the pertinent data is tranferred into the edited 
file. Inappropriate cards are removed. If the English units card 0001 is present 
in the edited file, the units will be adjusted from their default MKS values. The 
aircraft is assumed to be flying at a height of 10.0 units (meters or feet); the 
ata ment must be corrected after aircraft characteristics are transferred to the 

ile. 


Recover Drop Size Distribution From Database 


AGEDIT also contains access logic to FSCBG.MSL, the complete drop size 
distribution library from Skyler and Barry 1990. If the binary data file cannot 
be located, AGEDIT writes 


FSCBG.MSL Cannot Be Found Or Opened. 
Press Any Key To Continue ... 


after which AGEDIT returns to editing. 


Five drop size descriptions are shown on the screen at a time. Movement 
through the drop size distribution library is with the up and down cursor arrow 
keys. The selection of a specific drop size distribution is made by positioning 
the highlight line on the desired description and pressing the ENTER key. The 
menu option may be aborted with the F1 function key. 


When a drop size distribution is selected, the pertinent data is tranferred into the 
edited file. Inappropriate cards are removed. FSCBG.MSL contains complete 
drop size distributions that have been combined by pairs if more than 16 drop 
sizes were originally present. 

Recover Drop Size Distribution From SDC Savefile 


AGEDIT permits the importing of drop size distribution information from 
DROPSIZE (Teske 1990b) and SDC (Teske 1989d) with the message 


Enter SDC Savefile FILENAME.EXT 
If the file cannot be located, AGEDIT returns to editing. Otherwise, it opens 
the file, assumes it is in SDC file format, and reads and transfers drop size 
distribution data to the edited file. Only the first 16 drop sizes are transferred to 
the edited file. 
Compute Wet Bulb Temperature Depression 
AGEDIT contains a simplification of the WETBULB program (explained 
below) by computing the wet bulb temperature depression (for card 0065) with 
the following entries 

Enter Dry Bulb Temperature 

Enter Relative Humidity 
Temperature is deg C unless card 0001 is present in the edited file; then the 


temperature is deg F. Relative humidity is in percentage. Sea-level pressure is 
assumed. Results are computed and entered on card 0065 with a volatile 
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fraction of 0.875. The actual volative fraction must be corrected after the wet 
bulb temperature depression is transferred to the edited file. 


In all cases, a return to the edited file places the cursor in the home position at the 
top of the file. 
AGVIEW 


AGVIEW presents the plotted results in a PLT file with the command 
AGVIEW filename 

If a file name is not entered on the command line, AGVIEW asks 
Enter PLT File To View 

If the entered PLT file cannot be found, the message is 
filename.PLT Cannot Be Found Or Opened 

The actual plots as sequentially produced in AGPLOT are: 


1. First shown on the screen. Some of the titles and scales may appear scrunched 
because of the way the plot is presented. 


2. Offered for transference to a printer attached to the computer by responding Yes 
to the question 


Hardcopy Previous Plot To Printer (N/Y) 


AGVIEW works by reading specifications about the personal computer screen and 
the attached dot matrix printer from a default communications file called AGVIEW.DEF. If 
this file is not found, the message is 


AGVIEW.DEF Cannot Be Found Or Opened 


AGVIEW.DEF has a four-line default structure. For VGA screens and IBM/Epson 
printers it takes the form: 


11 0.749609 
9M bare pe So bp 3H ih, Be 
Pad sd bik WS 
025 ee 2 


Included on the distribution diskette is SAMPLE.PLT, a sample plot (Figure 8.2) 
that may be used to confirm screen and printer characteristics with the command 


AGVIEW SAMPLE.PLT 
The graphics information generated by AGPLOT and read by AGVIEW is written 


for display on a Tektronix 4025. On this terminal screen there are 640 pixels horizontally 
and 742 pixels vertically. Plots are constructed on their side. AGVIEW rotates the plots 
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ninety degrees for viewing on the screen, and keeps them on their side for plotting on the 
printer. All personal computer terminal screens have 640 pixels horizontally, but vary in 
vertical pixels: 

CGA 200 pixels vertically 

EGA 350 pixels vertically 

VGA 480 pixels vertically 

Details about the four lines in AGVIEW.DEF are as follows. 

Lene 1: 


This line tells AGVIEW what type of screen the computer has and its resolution. 
The three default values currently on the market are: 


CGA 2 0.311424 

EGA 9 0.546166 

VGA 11 0.749609 
Lines 2, 3 and 4: 

Every number (without quotation marks) and character (with quotation marks) is 
separated by a comma, and each line ends with a period. In this notation numbers are taken 
as ASCII values, to be invoked with the CHR$ function (the escape character is 
CHR$(27)), while characters are used as typed. 

Laine 2: 7d ioe pea ee a gd Ge 
IBM/Epson printers enter graphics mode with a command string 
CHR$(27);"3";CHRS$(N); 
to specify the vertical movement of the paper per line. The "3" invokes spacing of N/216 
inches, and in this application, the value of N is 24. Thus, the entire 742 pixels should be 
dot-matrix printed into a space equal to 
(742/8) * (24/216) = 10.3 inches 
Changing the value of N modifies the appearance of the resulting plot. 
The additional information on the second line 


CHR$(27);"O"; 


turns off the page skip-over-perforation at the end of sixty lines. This option is necessary 
because of the length of the plot. 


Pine 3: D7 ales 2852: 
IBM/Epson printers enter the number of dots per line with a command string 


CHR$(27);"L";CHR$(N1);CHR$(N2); 
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to specify the number of dots to be plotted per inch. The "L" invokes 120 dots per inch. 
The next two numbers give the total number of dots to be plotted on each line, using BASE 
256 arithmetic. Thus, since 640 dots are needed, a value of N1 = 128 and N2 = 2 gives 
128 + 2 (256) = 640 
Line 4: Dea 
IBM/Epson printers terminate graphics with the command string 
CHR$(27);"2"; 


A definition file for Hewlett-Packard DeskJet and LaserJet printers takes the 
following appearance: 


11 0.749609 
21 Ebe27 eatlOORe27. *r0A” HP 
27,"*b80W". 
DF ie BOGE: 

The HP after the period on Line 2 is essential. If the printer is one of these types, 
simply replace AGVIEW.DEF with this definition file and correct for the screen type. 
AGCTRL 

AGCTRL provides a shell program for managing the calculation from AGEDIT to 
AGDISP to AGPLOT to AGVIEW. The shell is invoked with the command 

AGCTRL filename 
If a filename is not entered on the command line, AGCTRL asks 
Enter Master Input File Name 
The process options available are the following: 
-1 Change Master Input File Name 
0  ExittoDOS 
1 AGEDIT To Edit Master Input File 
2 AGDISP To Integrate Equations 
3. AGPLOT To Construct Plots 
When AGCTRL asks 
Enter Option To Invoke 


the program always presents the next most logical process option in parenthesis. Pressing 
ENTER invokes this default. 


If AGDISP or AGPLOT generates errors during program operation, AGCTRL asks 
Error Exit From AGDISP/AGPLOT -- Press ENTER To Continue 
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to give enough time to read the screen and discover the mistake. 
- AGDISP is actually invoked with the command line 
AGDISP filename.INP 
while AGPLOT is actually invoked with the command line 
AGPLOT filename.LOG 


WETBULB 


The wetbulb temperature depression entry on card 0065 may be computed 
separately by invoking the program WETBULB. WETBULB asks for the following 
information, and then computes the temperature depression using the Carrier equation 
(Jennings and Lewis 1950): 

Dry Bulb Temperature (deg F or deg C) 

Height Above Sea Level (feet or meters) 

Relative Humidity (0 to 100) 

The U. S. Standard Atmosphere is assumed to apply, and the ASME Steam Tables 
relate water vapor pressure to temperature. The Carrier equation is solved by iteration. 


DROPSIZE 


The DROPSIZE program (Teske 1990b) permits recovery and manipulation (simple 
interpolation and extrapolation) of the complete drop size distribution database from Skyler 
and Barry 1990. 


SDC 


The SDC program (Teske 1989d) permits entry of drop size distribution data (such 
as from DROPSIZE) and combining of these data into reduced input data sets for entry into 
AGDISP. 


AGCTRL 


Figure 8-1 Program flowchart for operation of AGDISP and AGPLOT 
on personal computers. 
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9. Atmospheric Turbulence Level 


The second entry on card 0050 is the maximum value of background turbulence q?. 


This value influences the growth of the variance about the mean trajectories, and the 
standard deviation of the deposition. The crosswind inputs on card 0028 also produce a 
turbulence level that is added to the card 0050 level, since a logarithmic profile generates a 
turbulence level of: 


q2 = 0.845 





V(zr) i 
In ( 2/ Zo ) (39) 


where V(z,) is the velocity at height z, with surface roughness z,. Under locally neutral 


atmospheric conditions, the total wind velocity would be used to compute the total 
turbulence level by the above formula. If the atmosphere is calm, the turbulence level may 
be taken as 0.0 on card 0050. 


For atmospheric conditions other than neutral or calm, estimates of turbulence 
levels may be obtained by using stability categories (Pasquill and Smith 1983, Csanady 
1973, or using the Richardson number entry on card 0027). Table 9-1 presents these 
averaged atmospheric categories as a function of surface wind and temperature inversion, 
while Table 9-2 gives the range in values for turbulence as a function of surface wind and 
stability category. Selecting the simulation time-of-day and surface conditions permits the 
selection of the stability category, A through G. This category then provides the selection 
of turbulence values consistent with the stability. The Richardson number approach results 
in turbulence factors consistent with those shown in Figure 9-1. 


TABLE 9-1 


Stability Categories in Terms of Wind Speed, Insolation and State of Sky 
= (from Pasquill and Smith 1983) 











INSOLATION NIGHT 
Surface 
wind speed Thinly overcastor  < 3/8 
(m/sec) Strong Moderate _ Slight > 4/8 low cloud cloud 
<2 A A-B B G G 
2-3 A-B B C E F 
3-5 B B-C (6: D E 
5-6 C C-D D D D 
>6 G D D D D 


Notes: 1) For A-B take average values for A and B; etc. 


2) Strong insolation corresponds to sunny midday in midsummer, slight 
insolation to similar conditions in midwinter. Night refers to the period 
from one hour before sunset to one hour after dawn. The neutral category 
D should also be used, regardless of wind speed, for overcast conditions 
during day or night, and for any sky conditions preceding or following the 
night as defined above. 


TABLE 9-2 


Turbulent Intensities Near Ground Level 
(from Csanady 1973) 


Stability Category q? / V(z,)* 
extremely unstable 0365.--.1-2 
moderately unstable 0.145 - 0.365 

C-D-E near neutral 0.025 - 0.145 
moderately stable 0.015 - 0.025 
extremely stable 0:0" - 0.015 
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Appendix A: AGDISP Subroutines 


This section of the Appendix summaries the code comprising each subroutine in 
AGDISP. Figure A-1 presents a flowchart of AGDISP. Table A-1 highlights the common 
block variables used in the program. 


AGDISP 


AGBKG 


AGBOD 


AGBZD 


AGBZG 


is the mainline program computing material motion by the Lagrangian 
technique. Each data entry is checked and initialization is performed. Card 
order is maintained by a series of error flags L20, L23, L25, etc. Free format 
entry requires that all data cards contain the appropriate data; otherwise, an 
error will be detected. All input cards generate output to the terminal screen. 
Integration is then performed, after which the final ground deposition is 
computed, including the effects of evaporation. The generated plot files are 
used when plotting the trajectories and deposition with AGPLOT. 


evaluates the background velocity and turbulence level at a specified point. 
The decay rate, based on the time constant DTAU, is first computed. The 
simple evaporation model of Trayford and Welch 1977, and the drag law of 


Langmuir and Blodgett 1946 are used. The turbulence q? is then computed. 
The turbulence may come from superequilibrium or the fixed input value 
modified by other models in the flow. Alternately, the turbulence may be 
provided in a WAKE plot file. The material-turbulence correlations are then 
computed. The equations are a result of the assumption of a von Karman 
spectral distribution (Bilanin and Teske 1984). The constants XK1, XK2, 
XK3 take their limiting values whenever the time constants are within one 
percent of each other. The correlations < ux > and < uv > are then evaluated. 
Finally, the needed information at the specified point is saved in the vector 
DV. 


computes the rate of change in the wide-body area, as a function of distance 
behind the aircraft nose, and initializes the turbulent wake parameters when 
the wide-body effect (with cards 0080 and 0081) is invoked. A three-point, 
nonuniform derivative formulation is used, with correction at the end points 
(the nose and tail). 


evaluates the derivatives for the time-dependent solution of the Betz roll-up 
methodology. The variables are the radius of the vortex core and the 
circulation strength of the vortex. At initialization the radius is zero and the 
singularity is treated by taking the appropriate limit. 


initializes the Betz roll-up procedure using the user-inputted circulation 
distribution stored in common block BETZ. The first section of the 
subroutine computes the spatial derivative of the circulation along the wing. 
The circulation distribution is then reexamined to locate where the slope is 
maximum (these become positions where vortices start rolling up) and where 
the slope is minimum (these become the end positions of the vortex circulation 
patterns). The circulation pointers to start, maximum and end for up to four 
vortices are then set (more than four vortices invoke an error exit from 
AGDISP). The initialization is completed by evaluating the derivatives and 
initializing all of the parameters pointing in the AGDISP code to vortex center 
locations, strength, canopy effect and unrolled sheet effect. The code then 
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AGBZI 


AGCDS 


AGCRS 


AGDATI 
AGDAT2 
AGEQN 


AGGLQ 


AGINT 


AGMAT 


establishes a vortex-dependent time step DTV vector based on one percent of 
the roll-up time constant. 


integrates the Betz equations across the time step DELT entered from the 
AGDISP code. The first section of the subroutine establishes the time step 
and number of steps to integrate based on the computed values of DTV and 
the vortices not yet fully rolled up. Each vortex is treated separately by first 
predicting its new values of radius and circulation, then solving for the 
derivatives at these values and correcting the solution. The centroid of the 
vortex 1s computed based upon how much circulation has rolled into the 
vortex at each time step. The incremental movement of the vortex for each 
Betz time step is added to the current position of the vortex (this position may 
be influenced by other background features incorporated in AGDISP), and 
rolled-up vortices are flagged with zero sheet strength. Lastly, the unrolled 
sheet lengths are determined. 


writes the input data deck at the front of the log file. If card 0005 is present, 
AGCDS examines the casefile to locate the proper case, then overlay these 
cards with any cards in the input deck. The resulting run file will be read in 
AGDISP as the input runstream. 


computes the local value of turbulent energy associated with the discrete 
crosswind velocity profile shape entered with 0029 cards. A three-point, 
nonuniform derivative formulation is used to compute the velocity gradient. 
Locally neutral superequilibrium is assumed to compute the QV profile. 


processes half of the input cards. 
processes the other half of the input cards. 


monitors the integration of the equations. Initially, all positions are stored and 
initial derivatives determined. Integration then proceeds step by step to the 
maximum time. The WAKE plot file, Betz, canopy, propeller or jet engine, 
helicopter and wide body are updated where applicable. Solutions terminate 
at the deposition height. Since the vortices move under the influence of one 
another and their ground images, their positions are also adjusted. Any 
material that impacts the surface is flagged and its final size is computed. 
Termination is checked for evaporation, maximum time, and ground impact; 
and plot save is invoked. If termination exists, transfer is returned to 
AGDISP; otherwise another time step is taken. 


is a Gauss-Legendre integration routine used to determine the area under the 
discrete function between starting and ending points (Carnahan, Luther and 
Wilkes 1969). When M equals 1, the integrand is multiplied by y to generate 
the spatial moment of the discrete function. 


is a linear interpolator extracting the value of the discrete function at the 
desired position. 


fills the six-by-six matrix array for the unknowns < uu >, < vv >, < ww >, 


_ <uv >, < uw D, and < vw > by superequilibrium for a given value of q?. 


The linear equations are solved to determine the difference between q* and 
<uu>+<vw>+<ww>. 
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AGPAD 


AGSAV 


AGSRI 
AGSUP 


AGSVE 


AGVEL 


AGWAK 


computes the vortex circulation reduction resulting from interactions with a 
canopy. Essentially, the scrubbing of the vortex acts as a drag on the wake 
flow field. The drag translates into an effective vortex strength smaller than 
the noncanopy value. In this subroutine vortices are checked for whether they 
have penetrated the canopy; if so, an integration across the portion of the 
vortex interacting with the plant area is taken, a simple time integration is 
performed, and the vortex strength reduction factor is computed. The 
reduction factors FACR and FACL modify the vortex strengths, unrolled 
sheet strengths, propeller swirl and helicopter effects. 


writes the step integration results from AGDISP to the plot file (mean 
position, standard deviation, etc.) and the line printer (all variables for all 
nonimpacted material). With the proper flags set, this routine will also write 
to the plot file and print vortex center and powerplant locations. 


evaluates the turbulence correction for nonzero Richardson number. 


is the controlling routine for superequilibrium. Given the six spatial 
derivatives in the AGDISP code and the local scale length, the 
superequilibrium equations are iterated for the values of < uu >, < vv > and 
< ww >. The maximum gradient is used to normalize the solver, an 
approximate result is selected, and then an accurate solutions is found by 
stepping in q until a zero crossing is found by bisection. This result produces 
the values of turbulent energy unnormalized by scale length and the 
normalizing velocity gradient. 


computes the incremental background velocity from the unrolled Betz sheets. 
In this case the sheets are represented by a constant circulation sheet of 
vorticity, and the resulting flow is analytically determined. The singularity in 
the vicinity of the sheet is controlled by imposing a linear profile in z. As the 
Betz procedure continues, the sheet becomes shorter until it is finally rolled 


up. 


computes the mean velocity components U, V, and W at the position (x,y,z). 
Each vortex, its reflecting image across y = 0 and their images across z = 0, 
are used to compute the overall velocity increment. The standard potential 
vortex velocity field producing a velocity normal to the radius vector is broken 
into its (y,z) components and modified by FACR and FACL for the presence 
of a canopy and GDKYV for circulation decay. The unrolled sheet effect (with 
image sheets) is added, as are the helicopter, powerplant, wide body and 
mean crosswind modified by the canopy. Alternately, the WAKE plot file is 
quizzed to return the appropriate velocity components. 


processes calls for information in the WAKE plot file. The file is sequential 
binary, with a first record of the number of y and z mesh points, a second 
record of the values of the y mesh, and third record of the values of the z 
mesh. Following records consist of a time entry followed by all y data values 
for each z row for each variable. The second section recovers the WAKE plot 
file data for the desired time T. Interpolation between time steps is performed 
every one-tenth of the interval; otherwise the data arrays in common block 
WAKE are not updated. The file is rewound, the initial data skipped, and the 
pertinent time step bracketed. Linear interpolation is then performed, and the 
next time check is computed. When the end of the WAKE plot file is reached, 
the last time entry data is used for the duration of the AGDISP run. The third 
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section linearly interpolates the data array read from the WAKE plot file for 
the variable desired (V,W, q2) at the position (y,z). A warning message is 
written the first time spatial extrapolation beyond the WAKE plot file grid 
coordinate is invoked, even though extrapolation may continue indefinitely 
thereafter. 


AGVEL AGSUP AGWAK 
AGSVE AGMAT 


Figure A-1 AGDISP summary flowchart. 


Name 
AV 


CCW 
CEFE 
CHF 
CHG 
CHQ 
CHR 
CHS 
CHW 
CHZ 
CPUR 
CPXI 
CPXIS 


DBOD 
DCUT 
DEA 
DEAS 
DENF 


DIAM 
DIAMV 


Common 


SUPR 
WAKE 


HELI 
HELI 
BODY 
BODY 
BODY 
HELI 
HELI 
PROP 
PROP 
LOOP 


BODY 
EVAP 
EVAP 
LOOP 
EVAP 
BETZ 

EVAP 

DROP 


TABLE A-1: AGDISP Common Variable List 


Description 

vector: plant area fraction 

vector: superequilibrium values 

vector: wake plot file contents 

vector: Canopy parameter indicates amount of material removed 
cosine of crosswind angle 

canopy collection efficiency 

helicopter transition factor (hover to vortex pair) 
helicopter circulation factor (vortex pair) 

wide-body wake turbulence factor 

wide-body wake radius factor 

wide-body interpolated cross-sectional area 

helicopter downwash factor (hover) 

effective surface roughness 

powerplant wake turbulence factor 

powerplant virtual origin factor 

saved initial powerplant virtual origin factor 

cosine of terrain angle 

current diameter of wide body 

drop diameter below which evaporation ceases 
parameterized evaporation constants (also DEB and DEC) 
vector: parameterized evaporation constants (also DEBS and DECS) 
specific gravity 

vector: Betz circulation derivatives 

drop diameter 


vector: drop sizes 
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DIST 
DMASS 
DMCV 


DOVS 
DSYM 
DSYMS 
DSYP 
DSYPS 


DTAU 
DTEMP 
DTV 
DZBP 
EDNV 
EDOV 
FACL 
FACR 
FHEL 


FLOW 
FR 
G2PI 
G2PIS 
GDK 
GDKV 
GSAV 


INPT 

DROP 
EVAP 
VORT 
LOOP 


LOOP 


LOOP 
NORM 
NORM 
EVAP 
VORT 
INPT 
NORM 
NORM 
MEAN 


HELI 

INTG 
DROP 
INTG 


LOOP 
GAMD 
GAMD 


vertical distance to nominal height 

vector: mass fractions 

vector: drop volume ratio at surface impact 

vector: Betz roll up derivatives 

vector: saved initial Betz roll up derivatives 

vector: unrolled Betz sheet length left-of-centroid 

vector: saved initial unrolled Betz sheet length left-of-centroid 
vector: unrolled Betz sheet length right-of-centroid 

vector: saved initial unrolled Betz sheet length right-of-centroid 
time step 

drop relaxation time 

wet bulb temperature difference 

vector: Betz integration time step for each vortex 
incremental vertical distance from DIST to biplane wing 
vector: present drop diameter size 

vector: previous drop diameter size 

vector: plant canopy circulation reduction for left vortices 
vector: plant canopy circulation reduction for right vortices 
helicopter distance factor 

vector: plant area fraction integral for left vortices 

flow rate of spray system 

vector: plant area fraction integral for right vortices 

vector: circulation divided by two pi 

vector: saved initial circulation divided by two pi 

vortex decay constant 

vector: vortex decay factor for each vortex 


vector: Betz average sheet circulation per unit length 
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GSAVS 
GV 
HHEL 
HTPAD 
ICLOC 
ICMX 
ICV 
IDLOC 
IOUT 
ISW 
ISWC 
IWLOC 
IXLOC 
JHEL 
KVOR 
LBOD 
LEVAP 
LMCRS 


LPRP 
LQQSE 
LUNIT 
LV 

LX 


LOOP 
BETZ 
HELI 


SAVE 


SAVE 
SAVE 
SAVE 
SAVE 
SAVE 
SAVE 
HELI 

VORT 
BODY 
EVAP 


vector: saved initial Betz average sheet circulation per unit length 
vector: Betz circulation 

initial height of helicopter rotor plane above surface 
plant canopy maximum height 

flag to indicate saving of canopy information in plot file 
number of input cards 

vector: card numbers 

flag to indicate saving of drop size information in plot file 
flag to force output to screen 

vector: flag indicating ground encounter 

flag to sum number of drops hitting the ground 

flag to indicate saving of vertical velocities in plot file 
flag to indicate saving of all velocities in plot file 
dividing streamline flag 

flag indicating number of vortices not yet rolled up 
running pointer indicating position along wide body 
evaporation flag 

mean crosswind flag 

mean velocity flag 

flag indicating powerplant type (jet engine or propeller) 
turbulence flag 

English units flag 

vector: card flags 

error flag 

vector: end points location for all Betz vortices 

vector: starting points location for all Betz vortices 


vector: maximum points location for all Betz vortices 


N2FRO 
NBOD 


NCRS 
NDAT 
NDEF 
NDRP 
NENDF 
NEXTF 
NGAM 


NOUT 
NPAD 
NPLT 
NPRP 
NPRT 
NPTJ 


NPVX 
NPXX 
NSCR 
NSTEP 


NVOR 
NWPF 


NZ 


DROP 
BODY 
VORT 


DROP 
WAKE 
WAKE 
BETZ 


OUTP 


OUTP 
PROP 
OUTP 
OUTP 
OUTP 
OUTP 
OUTP 


SAVE 
OUTP 


WAKE 
WAKE 
WAKE 


"fro" flag 

number of discrete wide-body data points 
number of Betz vortices 

number of discrete crosswind velocity data points 
data file unit number 

casefile unit number 

number of drop sizes 

WAKE plot file record position counter 

WAKE plot file extrapolation flag 

number of discrete Betz circulation data points 
log file unit number 

screen Output unit number 

number of discrete plant area fraction data points 
plot output file unit number 

number of powerplant units 

printer output file unit number 

frequency of trajectory plot output file writes 
plot save flag for vertical velocity 

frequency of vortices/centerlines plot output file writes 
plot save flag for mean velocities 

temporary file unit number 

step number counter 

number of nozzles 

number of vortices 

WAKE plot file unit number 

number of horizontal points in WAKE plot file 
number of vertical points in WAKE plot file 


PGBP 
PSBP 
PV . 
QQBOD 
QQMX 
QV 
RBOD 


TMCV 


UO 
USK 


UZ 


BODY 
HELI 


PROP 
LOOP 
INPT 

BODY 


TFRO 
MEAN 
LOOP 


INPT 


NORM 


EVAP 
WAKE 


SUPR 
SUPR 


biplane wing circulation factor 

biplane wing semispan factor 

vector: vertical locations of crosswind velocity 
centerline turbulence level in wide-body wake 
constant background turbulence level 

vector: computed crosswind turbulence level 
current radius of wide body 

helicopter blade radius 

vortex core radius 

radius of propeller/jet flow field 

saved initial radius of propeller/jet flow field 
aircraft semispan/rotor radius 

vector: wide-body area 

sine of crosswind angle 

vector: rotation directions of propellers 
vector: vortex core radius 

vector: saved initial vortex core radius 

sine of terrain angle 

terrain slope angle 

final time value 

maximum integration time 

vector: time of impact with surface 

WAKE plot file time 

aircraft flight speed 

shear stress constant for crosswind evaluation 
vector: crosswind U velocity 


y derivative of U axial velocity 
z derivative of U axial velocity 
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VFRAC 
VPRP 
VPRPS 


XOS 
XOV 
XOVS 
XPRP 


YMSN 


DROP 
PROP 
LOOP 


SUPR 
SUPR 


HELI 


SUPR 
SUPR 


BODY 
BODY 
MEAN 
LOOP 
VORT 
LOOP 
PROP 
LOOP 


LOOP 


LOOP 
BODY 
LOOP 
HELI 

LOOP 


volatile fraction 

propeller tip speed 

saved initial propeller tip speed 
vector: crosswind V velocity 


y derivative of V horizontal velocity 
z derivative of V horizontal velocity 


helicopter downwash velocity at blade plane 


y derivative of W vertical velocity 
z derivative of W vertical velocity 


axial distance from trailing edge (shaft centerline) to wide body nose 
vector: wide body axial distance 

axial distance measured from nozzle 

saved initial axial distance measured from nozzle 
vector: Betz roll up values 

vector: saved initial Betz roll up values 

axial location of powerplants 

vector: saved initial nozzle locations 

vector: horizontal location of left vortices 

vector: saved initial horizontal location of left vortices 
vector: horizontal location of right vortices 

vector: saved initial horizontal location of right vortices 
horizontal location of wide-body centerline 

saved initial horizontal location of wide-body centerline 
horizontal location of helicopter centerline 

saved initial horizontal location of helicopter centerline 
minimum y value of mean trajectory 

maximum y value of mean trajectory 


minimum y value of spread trajectory 
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YMSX 
Yov 
YOVS 
YPRP 
YPRPS 


ZBALS 
ZBAR 
ZBARS 
ZBOD 
ZBODS 


ZHELS 
ZMSN 
ZMSX 
ZOPAD 
ZOVS 
ZPRP 


ZPRPS 
ZREF 


VORT 
LOOP 
PROP 
LOOP 


BETZ 
WAKE 


LOOP 


LOOP 
BODY 
LOOP 
HELI 

LOOP 


VORT 
LOOP 
PROP 
LOOP 
NORM 


WAKE 


maximum y value of spread trajectory 

vector: Betz roll up horizontal locations 

vector: saved initial Betz roll up horizontal locations 

vector: horizontal location of powerplant centerline 

vector: saved initial horizontal location of powerplant centerline 


vector: discrete stations for Betz circulation 
vector: horizontal points in WAKE plot file 


vector: vertical location of left vortices 

vector: saved initial vertical location of left vortices 
vector: vertical location of right vortices 

vector: saved initial vertical location of right vortices 
vertical location of wide-body centerline 

saved initial vertical location of wide-body centerline 
vertical location of helicopter centerline 

aed initial vertical location of helicopter centerline 
minimum z value of mean trajectory 

maximum z value of mean trajectory 

minimum z value of spread trajectory 

maximum z value of spread trajectory 

surface roughness height 

plant canopy effective surface roughness height 
vector: Betz roll up vertical locations 

vector: saved initial Betz roll up vertical locations 
vector: vertical location of powerplant centerline 
vector: saved initial vertical location of powerplant centerline 
apparent surface height 


vector: vertical locations of plant area fraction 
vector: vertical points in WAKE plot file 
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Appendix B: AGPLOT Subroutines 


This section of the Appendix summarizes the code comprising each subroutine in 
AGPLOT. Figure B-1 presents a flowchart of AGPLOT. Table B-1 highlights the 
common block variables used in the program. 


AGPLOT 


AGADD 


AGASC 


AGCAN 


AGCHK 


AGCHR 
AGCON 
AGCOV 


AGDEP 
AGDFT 
AGDRP 


AGEQD 


AGEQG 


is the mainline program for plotting the results from AGDISP. The program 
first reads the log file, recovers minimum and maximum values and identifies 
plot options. The user then selects a plot option, including plotting mean 
trajectories, trajectories including the turbulent standard deviation, trajectories 
of vortices or powerplant centerlines, ground deposition, equivalent density, 
droplet diameter or droplet velocity. More than one option may be invoked 
before program exit. 


computes a running total (as a function of time) of the deposition by adding 
contribution incrementally. 


computes minimum, maximum and increment scale values to include available 
data and present a visually pleasing scale division. The algorithm uses a base 
ten power law to decide on scale divisions. 


computes the incremental Gaussian deposition through a canopy. Each 
released material moves through the canopy, and its deposition effect is 
summed and contour plotted by the user. 

verifies the appropriateness of the scaling information requested by AGPLOT 
of the user. In all cases the scale must increase with a positive spacing 
between scale divisions and an integer number of scale divisions not to exceed 
ten. 

plots the character strings consistent with the option selected by the user. 
computes the continuous ground disposition pattern. 


computes the coefficient of variation as a function of lane separation between 
flight lines and overlaps deposition to generate the composite swath pattern. 


evaluates the incremental continuous deposition. 

computes the drift fraction. 

reads the selected plot file to recover the time history of a specified drop 
diameter undergoing evaporation, the canopy parameter or the mean velocities 
of a drop. 

computes the equivalent Gaussian distribution. The mean horizontal and 
vertical positions and standard deviations are summed over all material in the 
solution plane, and the figure of merit is computed. 


interprets the selected plot file to construct an equivalent Gaussian profile 
distribution at every time step. Two passes are made through the plot file. 
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The first pass computes the equivalent Gaussian, figure of merit and vertical 
velocity criterion at every time in the plot file, and displays to the screen the 
points of maximum and minimum figure of merit, material impact with the 
surface, and equivalent Gaussian impact with the surface. From this 
information the maximum figure of merit is determined and the program 
locates that time in the second pass to create the equivalent plot. 


AGEQP sets up and plots the equivalent Gaussian and material isopleths. 
AGIMP computes the velocity normal to the collection device, the impaction 
. parameter, impaction efficiency and incremental material deposition in a time 
step. 
AGNUM _ converts a real number into a character string. 
AGOPN opens a selected BIN plot file. 


AGPLT is the plotting subroutine that uses graphics subroutine calls to perform the 
actual plotting. The screen is first cleared, then axes are established and 
plotted and a wing schematic is drawn. The plot file is read so that nine 
points are plotted at a time, drawing the mean trajectories, standard deviation 
paths, vortices and powerplant centerline paths, tag markers at a desired 
interval, array profiles (including ground deposition), and equivalent 
Gaussian profiles. 


AGRPF rereads the front of the log file to recover input card data for crosswind 
velocity or plant area profiles. 


AGSET initializes scale sizes for plotting options. 


AGSPD uses the stored plot file variables to determine the plotting coordinates of the 
standard deviation, normal to the instantaneous direction. 


AGSRF oversees computation of Gaussian deposition on-a specified collection device 
at a specified location (y,z), with a specified unit normal. 


AGSRT determines the scaling to recover the desired deposition units. 
AGSUM _seintegrates the deposition along the surface. 


AGVMD computes the volume and number median diameters for the selected ground 
deposition. 


AGVRE computes the flux of material through a specified vertical surface. 


Graphics calls on the Data General invoke GKS subroutine calls. They include the 
following: 


GKS_CLEAR clears the graphics screen. 
GKS_CLOSE ends the graphic session. 
GKS_DRAW connects the data points with straight lines after setting up the axes. 
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GKS_GPL 
GKS_OPEN 
GKS_Q_WS_WIN 
GSK_QPRINT 
PLDG 


plots the data points. 

initializes the graphics session. 

returns the pixel size of the plotting device. 
sends hardcopy to the output device. 


plots pixel data. 


Screen calls on the Data General invoke the following subroutine calls: 


FREAD 

IREAD 
SCREEN_INIT 
SCROLL 
SREAD 
SSCREEN 
SWRITE 


recovers real number information from the user. 
recovers integer number information from the user. 
initializes the screen. 

perform a scroll of the data screen. 

recovers character information from the user. 
rapidly replaces screen contents. 


writes information to the screen. 


Graphics calls on VAX and on IBM PC/XT/ATs are invoked with Tektronix 4025 
calls. These calls write plot data to the screen or a data file for later viewing. The 
subroutines invoked during the running of AGPLOT include the following: 


PLHC 
PLLN 
PLNT 
PLPL 
PLSC 
PLST 
PUTT, 


completes the plot. 

sets the line type. 

sets up scale notation and labels. 

plots the data points (connected by straight lines). 
initialized the graphics area and scale increments. 
writes a character string. 


establishes the default graphics parameters. 


AGPLOT 
AGCAN AGDFT AGEQG AGSRF 


_ Figure B-1_ AGPLOT summary flowchart. 


B-4 


Name 
AV 
CBAR 
CNUM 
CVEC 


DELT 
DIAM 
DIAMT 
DIAMV 
DIST 
DMASS 
DY 
DZBP 
FLOWR 
FNAME 
FOM 
GGRD 
GTEM 
GVEC 
HTPAD 
ICOL 
ICV 
IOPT 
ISW 


Common 
PVEC 
EQGN 
NUMB 
PVEC 
EQGN 
PLOT 
SURF 
SURF 
DROP 
PLOT 
DROP 
CONT 
PLOT 
DROP 


EQGN 
STOP 
PVEC 
AVEC 
PLOT 
SURF 
EQGN 
PLOT 
PVEC 


TABLE B-1: AGPLOT Common Variable List 


Description 

vector: solution values for plotting 

equivalent Gaussian magnitude 

displayed number for several plot options 

vector: VMD/NMD/canopy intermediate values 
vector: drop diameter factors for equivalent Gaussian 
tag time increment 

drop size 

collector significant dimension 

vector: drop diameter 

initial vertical distance from surface to aircraft wing 
vector: mass fraction 

y horizontal scale increment 

incremental vertical distance from DIST to biplane wing 
spray system flow rate 

master file name 

equivalent Gaussian figure of merit 

vector: ground and canopy deposition values 
vector: temporary value array for vertical flux 
vector: multiple array vector 

canopy height - 

flag indicating collector type 

vector: surface impact flag 

plot option number 


vector: surface impact flag 


ITYPS 
IVAR 
IVEC 
JNAME 
KGRD 
KTYP 
KVOL 
L2FRO 
LAPLT 
LCLOC 
LDLOC 
LDROP 
LEPLT 
LGRID 
LMPLT 
LNVOL 
LPRP 
LSCAL 
LSPLT 
LTPLT 
LVEL 
LVOR 
LVPLT 
LWLOC 
LXLOC 
N2FRO 


AVEC 
NUMB 
PLOT 
PLOT 


STOP 
STOP 
STOP 
DROP 
PLOT 
FLAG 
FLAG 
DROP 
PLOT 
DROP 
PLOT 
DROP 
PLOT 
DROP 
PLOT 
PLOT 
PLOT 
PLOT 
PLOT 
FLAG 
FLAG 
DROP 


deposition scale flag 

flag indicating CNUM type 

nozzle number flag 

drop velocity flag 

number of characters in FRAME 

vector: number of entries in YGRD and GGRD 
vector: flag indicating ITYP for GGRD 

vector: flag indicating volatile/nonvolatile for GGRD 
flag indicating selection of "to" or "fro" direction 
array plot flag 

flag indicating canopy information in plot file 

flag indicating drop size information in plot file 

flag indicating selection of drop size 

equivalent Gaussian distribution plot flag 

flag indicating selection of grid autoscale 

mean trajectory plot flag 

flag indicating selection of volatile or nonvolatile contribution 
flag for powerplant centerlines to plot 

flag indicating selection of deposition scale units 
standard deviation trajectory plot flag 

tag plot flag 

number of variables present on plot file 

flag for vortex centroids to plot 

vortices/powerplant trajectory plot flag 

flag indicating saving of vertical velocities on plot file 
flag indicating saving of all velocities in plot file 


"fro" direction flag 


NDMN 
NDMX 
NDRP 
NINU 


NPLT 
NOUT 
NSAV 


NVEC 


PV 


SPV 


TA 
TEMND 


YCB 
YCV 
YGRD 


DROP 
DROP 
DROP 
OUTP 
OUTP 
OUTP 
OUTP 
OUTP 
DROP 
AVEC 
CONT 
PVEC 

PLOT 

EQGN 
PLOT 

CONT 
DROP 
PSPD 

PLOT 

DROP 
PVEC 

PVEC 

SURF 
SURF 
EQGN 
EQGN 
STOP 


minimum drop size file indicator 
maximum drop size file indicator 
number of drop sizes 

terminal input unit number 

log file unit number 

plot input file unit number 

screen output unit number 

graphics output unit number 

number of nozzles 

number of points in array plot vectors 
drop size increment pointer 

vector: powerplant locations for plotting 
aircraft semispan/rotor radius 

vector: standard deviation for equivalent Gaussian 
terrain slope angle 

scale factor 

vector: maximum simulation times 
tangent of terrain angle 

aircraft flight speed 

volatile fraction 

vector: left vortex locations 

vector: right vortex locations 

collector axial surface normal 

collector horizontal location 

equivalent Gaussian horizontal location 
vector: horizontal centerline locations for equivalent Gaussian 


vector: ground and canopy deposition horizontal scales 


B-7 


YMSNV 
YMSXV 


YSB 


VEC 
LE 
ZCB 


ZMSNV 
ZMSXV 
ZN 
ZREF 
ZSB 
ZTEM 
ZVEC 


PLOT 
PLOT 
DROP 
DROP 
CONT 
DROP 
DROP 
SURF 
EQGN 
PVEC 
AVEC 
SURF 
EQGN 
EQGN 
PLOT 
PLOT 
DROP 
DROP 
DROP 
DROP 
SURF 
PLOT 
EQGN 
PVEC 
AVEC 


horizontal axis plot scale maximum 

horizontal axis plot scale minimum 

vector: minimum horizontal value for mean trajectories 
vector: maximum horizontal value for mean trajectories 
YMIN for continuous deposition 

vector: minimum horizontal value for spread trajectories 
vector: maximum horizontal value for spread trajectories 
collector horizontal surface normal 

equivalent Gaussian horizontal standard deviation 
vector: temporary horizontal array for vertical flux 
vector: horizontal array vector 

collector vertical location 

equivalent Gaussian vertical location 

vector: vertical centerline locations for equivalent Gaussian 
vertical axis plot scale maximum 

vertical axis plot scale minimum 

vector: minimum vertical value for mean trajectories 
vector: maximum vertical value for mean trajectories 
vector: minimum vertical value for spread trajectories 
vector: maximum vertical value for spread trajectories 
collector vertical surface normal 

surface reference height 

equivalent Gaussian vertical standard deviation 

vector: temporary vertical array for vertical flux 


vector: vertical array vector 
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